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JET FORMATION IN BLACK HOLE ACCRETION SYSTEMS I: THEORETICAL UNIFICATION MODEL 

Jonathan C. McKinney' 
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ABSTRACT 

Two types of relativistic jets are suggested to form near accreting black holes: a potentially ultrarelativistic 
Poynting-dominatedjet and aPoynting-baryonjet. One source of jet matter is electron-positron pair production, 
which is driven by neutrino annihilation in GRBs and photon annihilation in AGN and x-ray binaries. GRB 
Poynting-dominated jets are also loaded by electron-proton pairs by the collisional cascade of Fick-diffused 
free neutrons. We show that, for the collapsar model, the neutrino-driven enthalpy flux (classic fireball model) 
is probably dominated by the Blandford-Znajek energy flux, which predicts a jet Lorentz factor of F ~ 100- 
1000. We show that radiatively inefficient AGN, such as M87, are synchrotron-cooling limited to F ~ 2 - 10. 
Radiatively efficient x-ray binaries, such as GRS1915H-105, are Compton-drag limited to F < 2, but the jet may 
be destroyed by Compton drag. However, the Poynting-baryon jet is a collimated outflow with F ^ 1-3. The 
jet from radiatively efficient systems, such as microquasar GRS1915H-105, may instead be a Poynting-baryon 
jet that is only relativistic when the disk is geometrically thick. In a companion paper, general relativistic 
hydromagnetic simulations of black hole accretion with pair creation are used to simulate jet formation in 
GRBs, AGN, and x-ray binaries. 

Subject headings: accretion disks, black hole physics, galaxies: jets, gamma rays: bursts. X-rays : bursts, 
supernovae: general, neutrinos 



1. INTRODUCTION 

Jets are a common outcome of accretion, yet the observed 
jet properties, such as collimation and speed, are not uniform 
between systems. This is despite the fact that the basic physics 
(general relativistic magnetohydrodynamics (GRMHD)) to 
describe such systems is black hole mass-invariant. Thus, it is 
worth-while to determine the unifying, or minimum number 
of, pieces of physics that would explain most of the features 
of gamma-ray bursts (GR Bs), x-ray binaries, and active galac- 
tic nuclei (AGN) dGhisellini & Celotti200i lGhiseninil2003t 
I iMeier 2003}). To understand jet formation requires at least 
explaining the origin of the energy, composition, collimation, 
and Lorentz factor The goal of this paper, and the compan- 
ion numerical models paper (McKinnev 2005b), is to explain 
these for GRBs, AGN, and x-ray binaries. 

Primarily we discuss two types of jets: Poynting-dominated 
jets typically dominated in energy flux by Poynting flux and 
dominated in mass by electron-positron pairs for AGN and 
x-ray binaries, while dominated in mass by electron-proton 
pairs for GRBs ; and Poynting-baryon jets with about equal 
Poynting flux and rest-mass flux and dominated in mass by 
baryons. The latter are sometimes referred to as coronal out- 
flows due t o their or igin. Generically this model is simi- 
lar to, e.g., ISol et alJ ( 119891) . while here the motivation is 
based upon the results of recent GRMHD numerical mod- 
els. This two-component jet model is one key to under- 
standing the diversity of jet observations. The Poynting- 
dominated jet is likely powered by the Blandford-Znajek ef- 
fect, while the Poynting-baryon jet is likely powered by both 
Blandford-Znajek power and the release of disk gravitational 
binding energy (McKinnev 2005a). Collimation of the polar 
Poynting-dominatedjet may be due transfield balance against 
the broader Poynting-baryon jet or by self-collimating hoop 
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Stresses. 

Among all the black hole accretion systems, it appears that 
the least unifiable is the observed emission. While the ra- 
diative physics is not black hole mass invariant, the observed 
differences suggest that the environment likely plays a signifi- 
cant role in the emission. For example, while both blazars and 
GRBs exhibit non-thermal emission, long-duration GRBs are 
harder with highe r luminosity, while blazars a re softer with 
higher luminosity (Ghirlanda et al."2004','200f). Also, GRBs 
lead to apparently most of the energy in 7-rays and less than 
10% to the sub-7-ray afterglow (Piran 2005). On the contrary, 
blazars apparently release only 10% in 7-rays and the rest is 
produced in the radio lobe (Ghisellini & Celottill2002^ . De- 
spite the difficulties in understanding the emission processes 
in some jet systems, the jet itself is probably produced by a 
universal process. 

The disk and jet radiative physics are keys to understand- 
ing the evolution of the jet and why different systems have 
different terminal Lorentz factors. Through radiative anni- 
hilation of photons in AGN and x-ray binary systems, the 
radiative physics may illuminate the origin of jet composi- 
tion by determining the electron-positron mass-loading the 
Poynting-lepton jet, and so the Lorentz factor of the jet. 
For GRBs, the radiative annihilation of neutrinos and the ef- 
fect of Fick diffusion by free neutrons from the corona into 
the jet dLe Vinson & Eichleri 120031) may give an understand- 
ing of the Lorentz factor of the jet and the origin of baryon- 
contamination. 

The rest of this section briefly reviews the types of black 
hole accretion systems and discusses jets in each. At the end 
is an outline of the paper 

1.1. GRBs 

Neutron stars and black holes are associated with the 
most violent of post-Big Bang events: supernovae and 
some gamma-ray bursts (GRBs) and probably some x-ray 
flashes (XRFs) (for a general revi ew see .Wooslevi 119931 
IWheeler. Yi. Hoflich. & Wandl2000h . Observations of a su- 
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pernova light curve (SN2003dh) in the afterglow of GRB 
030329 suggest that at least some long-duratio n GRBs are 
probably associated with core-collapse events ( Stanek et alJ 
12003; Kawabata et al. 2003; Uemura et al. 2003; Hiorth et aQ 

Neutrino processes and magnetic fields are both important 
to understand core-collapse. In unraveling the mechanism 
by which core-collapse supernovae explode, the implemen- 
tation of accurate neutrino transport has been realized to be 
critical to whether a supernova is produced in simulations 
dMesser et al. and col laborators 1998). This has thus far 
been interpreted to imply that highly accurate neutrino trans- 
port physics is required, but this could also mean additional 
physics, such as a magnetic field, could play a significant role. 
Indeed, all core-collapse events may be powered by MHD 
proce sses rather than neutrino processes ( Leblanc & Wilson 
1970t ISv mbalistv 1984; Wooslev & Weaver 198^ 
^ncan & Thompso n .1992; .Khokhlov et aL .19991 
Akivama. Wheeler. Meier. & Lichtenstad3 120031) . 
Core-collapse involves shearing subject to the 
Balbus-Hawley instability as in accretion disks 
((Akivama, Wheeler, Meier, & Lichtenstadt 2003). All 
core-collapse explosions are significantly polarised, 
asymmetric, and often bi-polar indicating a strong 
role of rotation and a magnetic field (see, e.g., 
Wang & Wheeler 1996; Wheeler. Yi. Hoflich. & Wang 
2000; Wang. Howell. Hoflich. & Wheeler 2001; Wang et al. 
2002; Wang. Baade. Hoflich. & Wheeler. .2003.. and 'refer- 
ences therein). Possible evidence for a ma gnetic dominated 
outflow has been found in GRB 021206 (Coburn & Boggs' 
1^03), marginally cons istent with a magnetic outflo w directly 
from the inner engine ( iLvutikov. Pariev. & Blandfordl2003l) . 
although these observations remain controversial. 

Black hole accretion is the key source of energy for many 
GRB models. Collapsar type models suggest that a black 
hole forms during the core-collapse of some relatively rapidly 
rotating massive stars. The typical radius of the accretion 
disk hkely determines the duration of long -dura tion GRBs 
JWooslevI fT993l iPaczvnskil \l99§L iMacFadven & WooslevI 
l1999h . An accretion disk is also formed as a result of a neu- 
tron star or black hole collisions with another stellar object 
JNaravan et al. 1992, 2001). 

GRBs are believed to be the result of an ultrarelativistic jet. 
Indirect observational evidence of relativistic motion is sug- 
gested by afterglow achromatic light breaks and the "com- 
pactness problem" suggests GRB material must be ultrarel- 
ativistic with Lorentz factor T > 100 to emit the observed 
nonthermal 7-rays (see, e.g., iPiranii2005<) . Direct observa- 
tional evidence for re lativistic motion comes from radio scin- 
tillation of the ISM ('Go odmanll 199 7*) and measurements of 
the after glow emitting region from GRB030329 (Tavlor et a^ 
I2004albl) . 

Typical GRB jet models invoke either a hot neutrino-driven 
jet or a cold Poynting flux-dominated jet, while both allow 
for comparable amounts of the accretion energy to power the 
jet (Popham et al. 1999). A neutrino-driven jet derives its en- 
ergy from neutrino annihilation from gravitational energy and 
the jet is thermally accelerated. However, strong outflows 
can be magnetically d riven (Bisnovatvi-Kogan & Ruzmaiki^ 
Il976t 'Lovelace" 1976^. 'Br andford.1976,) . In particular, black 
hole rotational energy can be extracted as a Poynting outflow 
llBiandford & Znaiekl 19771) . 

1.2. X-ray Binaries 



Long after their formation, neutron stars and black 
holes often continue to produce outflows and jets 
jMirabel & Rodriguez! fT99l These include x-ray binaries 
(for a review see lLewiiieralJl 995; McClintock & Remillardj 
2003), neutron star as pulsars (for a review see LorimM 
2001 on ms pulsars and Thorsett & Chakrabartv 19991 
on radio pulsars) and soft -gamma ray repeaters (SGRs ) 
jThompson & DuncanlfT99^ Il996t iKouveUotou etai]rT99l . 
In the case of x-ray binaries, the companion star's solar- 
wind or Roche-lobe forms an accretion disk. Many x-ray 
binaries in their hard/low state (and radio-loud AGN) show 
a corre lation between the x -ray luminosity and radio lumi- 
nosity (iMerloni et alJl20()3l) . which is consistent with radio 
synchrotron emission from a jet and x-ray emission from a 
geometrically thick, optically thin, Comptonizing disk. 

Some black hole x-ray binaries have jets (Mirabel et alJ 
[^92; Fender 2003), such as G RS 1915-H105 with apparently 
superluminal motion (F ^ 3) ('Mirabel & Rodriguez 19941 
IMirabel & Rodri guez jl999t Ipe nder & Belloni 2004), but may 
have F 1 .5 ( Kaiser et al. 2004). Synchrotron radiation from 
the jet suggests the presence of a magnetized accretion disk. 
Observations of a broad, shifted, and asymmetric iron line 
from GRS 1915H-105 is possibl e evid ence for a relativistic 
accretion disk ( Martocchia et a QEool, although this feature 
could be produced by a jet component. 

The standard paradigm is that relativistic jets from x-ray 
binaries are probabl y produced by the Blandford-Znajek ef- 
fect. However. .Gierlinski & Dona ( l2004) suggest that at least 
some black holes, such as GRS 1915H-105, have slowly rotat- 
ing black holes. If this is correct, then another mechanism is 
required to produce jets. Indeed, jets or outflows are produced 
from systems containing NSs, young stellar objects, super- 
soft x-ray white dwarfs, symbiotic white dwarfs, and even UV 
line-driven outflows from massive O stars. Indeed, a baryon- 
loaded coronal outflow with F^ 1.5-3canbe produced from 
a black hole accretion disk and not require a rapidly rotat- 
ing black hole ( Mc Kinnev & Gamm ie 2004). Nonrelativis- 
tic outflo ws were found eve n in viscous hydrodynami c sim- 
ulations dStone et al. I Il999[ [Tg umenshchev & Abram owicz I 
fT999l 1200(1 iMcKinnev & Gammiei 12002) . Such baryon- 
loaded outflows or jets are sufficient to explain most known 
x-ray binaries without invoking rapidly rotating black holes, 
and thus unifies such mildly relativistic jets in neutron star and 
black hole x-ray binaries. 

1.3. AGN 

Active galactic nuclei (AGN) have long been believed 
to be powe red by accretion onto supermassive black holes 
(Zerdovichi ll964l; ISatoeteil 11964 . Observations of MCG 
6-30-15 show an iron line feature consistent with emission 
from a rel ativistic disk with v/c ^ 0.2 (Tanakaet al. 199^ 
iFabian eralJl2002l) . although the lack of a temporal correla- 
tion between the continuum emissio n and iro n-line emission 
may suggest it is a jet-related feature ( lElvisl20 00). 

AGN are observed to have jets with F < 10 
( Urry & Padovani 1995; Biretta et al. 199 ^, even 
F ~ 30 (Begelmanet alJ 11994 GhiselH ni & Celotti 
l2(Kyll IW ad et al. 200 1 |), whi l e som e observatio ns im- 
ply F < 2 00 (Ghiselli nTgtai] 119931: iKrawczvns ki et alJ 
l2002HKonopelko et a l. 2003V Some radio-quiet AGN show 
evidence of weak jets (Ghisellini et al. 2004), which could be 
explained as a coronal outflow ( McKinney & Gammie 2004 
and not require a rapidly rotating black hole. Observations 
imply the existence of a two-component jet structure with a 
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Poynting jet core and a dissipative surround ing co mponent 
fchisellini & Madau 1996; Ghisellini et alj 120051) . The 
energy structure of the jet and wind are important in under- 
standing the feedback effect that controls size of the black 
hole and may determine the M -a relation dSpringel et alJ 
l2004^.ni Matteo et al. 2005) . 

1.4. Outline of Paper 

§ 121 summarizes the proposed unified model to explain jet 
formation in all black hole accretion systems. 

§|3]discusses why ideal MHD must break down in magneto- 
spheres and why the Goldreich-Julian charge density is never 
reached. A preliminary model is derived that describes the 
pair-loading and baryon-loading of the Poynting-dominated 
jet. 

§ 13 determines the Lorentz factor of Poynting-dominated 
jets. The GRB jet Lorentz factor is shown to be based upon 
electron-positron pair and baryon loading. We show that 
GRBs likely have electron-proton jets with Lorentz factor at 
large distances of 100 < Too < lOl 

Based upon pair creation rates for AGN and x-ray binaries, 
we show that relativity low radiatively efficient AGN, such as 
M87, have electron-positron jets with 2 < Too ^ 10. Radia- 
tively efficient systems, such as microquasar GRS1915H-105, 
likely do not have Poynting-dominated lepton jets but rather 
the observed jets are a relativistic coronal outflow from the 
inner-disk. 

§ |5] discusses Poynting-baryon jets and how they can ex- 
plain various observational features of jets in AGN and x-ray 
binaries. 

§|6l summarizes the key results and fits from GRMHD nu- 
merical models (McKinnev 2005b) used in this paper 
§0discusses the results and their possible implications. 
§|8]summarizes the key points. 

Appendix |a1 discusses breakdown of ideal MHD by 
electron-positron pair creation by radiative annihilation and 
electron-proton pair creati on by ambipolar and Fick diffusion. 
See also the discussion in iMcKinnevI (12005b ). Appendix 151 
gives a succinct summary of conserved flow quantities in 
GRMHD used in section^] Appendixiclgives a derivation for 
the lab frame stationary GRMHD forces along and perpendic- 
ular to the flow (field) line in the lab frame. This elucidates 
the origin of acceleration and collimation. Appendix|D]gives 
the formulae for Comptonization and pair annihilation used in 
section|3 

2. GRMHD PAIR INJECTION MODEL OF JET FORMATION 

The jet energy, composition, collimation, and Lorentz fac- 
tor are likely determined in a similar way for all black 
hole accretion systems. The particle acceleration mecha- 
ni sm and particle composition of the jet remained unexplained 
in iMcKinnev & Gammid (120041) . However, if field lines tie 
the black hole to large distances, then the source of mat- 
ter is likely pair creation since the amount of matter that 
diffuses across field lines is much smaller ('Phinne>ill2Sl 
iLevinson & Eichler 1993; Punslv 2001). Thus, the Poynting- 
dominated jet composition is electron-positron pair domi- 
nated in AGN and x-ray binaries. 

However, in GRB systems, free neutrons lead to baryon 
contamination due to Fick diffusion across the field lines and 
subsequent rapid collisio nally-induced avalanche decay to an 
electron-proton plasma ( Levinson & Eichler 2003). The pair 
annihilation rates are much faster than the dynamical time, 
and due to the temperature decrease, the electron-positron pair 
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Fig. 1. — Schematic of pair-production model and subsequent magnetic 
fireball formation for GRB disks. Fireball is extremely optically thick. Below 
a stagnation surface, pairs are accreted by the black hole and so do not load 
the jet. Here«g = GM/f-. 

rest-mass exponentially drops beyond the fireball formation 
near the black hole. Thus, the GRB jet composition is likely 
dominated by electron-proton pairs. 

GRMHD numerical models confirmed that accretion of a 
thick disk with height (H) to radius {R) ratio of ////?> 0.1 
with a homogeneous poloidal orientation self-consistently 
creates large scale fields that tie the black hole to larg e dis- 
tances (Mc Kinnev & G ammie 2004;"Hirose et al." 2004). Ac- 
cretion of an irregular field loads the jet with baryons and 
lowers the speed of the jet. However, the existence of a 
mostly uniform field threading the disk arises naturally dur- 
ing core-collapse supernovae and NS-BH collision debris 
disks. In AGN and solar-wind capture x-ray binary systems, 
the accreted field is probably uniform (Naravanet al. 2003f 
Uzdenskv & Spruit 2005). Roche-lobe overflow x-ray bina- 
ries, however, might accrete a quite irregular field geometry. 
The field geometry that arrives at the black hole, after travel- 
ling from the source of material (molecular torus, star(s), etc.) 
to the black hole horizon, likely depends sensitively on the 
reconnection physics. 

The reason why each system has some observed Lorentz 
factor has not been well-understood. One key idea of this 
paper is that the terminal Lorentz factor is determined by 
the toroidal magnetic energy per unit pair mass density en- 
ergy near the location where pairs can escape to infinity (be- 
yond the so-called "stagnation surface"). Put another way, the 
Lorentz factor is determined by the energy flux per unit rest 
mass flux for the rest-mass flux in pairs beyond the stagna- 
tion surface. For GRBs, neutron diffusion is crucial to ex- 
plain (and limit) the Lorentz factor For AGN and x-ray bi- 
naries, since a negligible number of baryons cross the field 
lines, pair-loading is crucial to determine the Lorentz factor 
of the Poynting-dominated jet since this determines the rest- 
mass flux or density. 

Figure ^ shows the basic picture for GRB systems, while 
figure 121 shows the basic picture for AGN and x-ray binary 
systems. An accreting, spinning black hole creates a mag- 
netically dominated funnel region around the polar axis. The 
rotating black hole drives a Poynting flux into the funnel re- 
gion, where the Poynting flux is associated with the coiling 
of poloidal magnetic field lines into toroidal magnetic field 
lines. The accretion disk emits neutrinos in a GRB model (7- 
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Fig. 2. — Schematic of pair-production model and subsequent shock- 
lieating and emission. AGN jet is optically thin and emits nonthermal and 
thermal synchrotron, while x-ray binary jet can be marginally optically thick 
and emit via self-absorbed synchrotron and by severe Compton drag. Severe 
Compton drag can lead to destruction of the Poynting-dominated jet. 



ray and many soft photons for AGN and x-ray binaries) that 
annihilate and pair-load the funnel region within some "injec- 
tion region." For GRB systems, neutrons Fick-diffuse across 
the field lines and collisionally decay into an electron-proton 
plasma. 

Many pairs (any type) are swallowed by the black hole, but 
some escape if beyond some "stagnation surface," where the 
time-averaged poloidal velocity is zero and positive beyond. 
Pairs beyond the stagnation surface are then accelerated by 
the Poynting flux in a self-consistently generated collimated 
outflow. In the electromagnetic (EM) jet, the acceleration pro- 
cess corresponds to a gradual uncoiling of the magnetic field 
and a release of the stored magnetic energy that originated 
from the spin energy of the black hole. 

One key result of this paper is that the release of magnetic 
energy need not be gradual once the toroidal field dominates 
the poloidal field, in which case pinch (and perhaps kink) in- 
stabilities can occur and lead to a nonlinear coupling (e.g. a 
shock) that conv erts Poynting flux into enthalpy flux ( Eichler 
ll993tlBegelmanH99 8). In the proposed GRB model, this con- 
version reaches equipartition and the jet becomes a "magnetic 
fireball," where the toroidal field instabilities drive large vari- 
ations in the jet Lorentz factor and jet luminosity. 

In AGN systems, nonthermal synchrotron from shock- 
accelerated electrons and some thermal synchrotron emission 
releases the shock energy until the synchrotron cooling times 
are longer than the jet propagation time. For AGN, jet ac- 
celeration is negligible beyond the extended shock zone, as 
suggested for blazars beyond the "blazar zone" dSikora et alJ 
I2005J . In x-ray binary systems, the shock is not as hot and 
also unlike in the AGN (at least those like M87) case the jet 
can be optically thick. Thus these x-ray binary systems self- 
absorbed synchrotron emit if they survive Compton drag. 

For all these systems, at large radii patches of energy flux 
and variations in the Lorentz factor develop due to toroidal 
instabilities. These patches in the jet could drive internal 
shocks and at large radii they drive external shocks with the 
surrounding medium. The EM jet is also surrounded by a 
mildly relativistic matter coronal outflow/jet/wind, which is 
a material extension of the corona surrounding the disk. This 
Poynting-baryon, coronal outflow coUimates the outer edge of 



the Poynting-dominated jet, which otherwise internally col- 
limates by hoop stresses. The luminosity of the Poynting- 
baryon jet is determined, like the Poynting-dominated jet, by 
the mass accretion rate, disk thickness, and black hole spin. 

This model is studied analytically in this paper, while in a 
companion paper we study this model numerically using ax- 
isymmetric, nonradiative, GRMHD simulations to study the 
self-consistent process of jet formation from black hole ac- 
cretion systems (McKinnev 2Q05b). Th ose simulations ex- 
tend the work of McKinnev & Gammie (2004) by including 
pair creation (and an effective neutron diffusion for GRB-type 
systems) to self-consistently treat the creation of jet matter, 
investigating a larger dynamic range in radius, and presenting 
a more detailed analysis of the Poynting-dominated jet struc- 
ture. 

Unless explicitly stated, the units in this paper have GM = 
c = 1, which sets the scale of length {rg = GM/c^) and time 
(tg = GM/c^). The mass scale is determined by setting the 
(model-dependent) observed (or inferred for GRB-type sys- 
tems) mass accretion rate {M[g i"' ]) equal to the accretion rate 
through the black hole horizon as measured in a simulation. 
So the mass is scaled by the mass accretion rate at the horizon, 
such that pojisk = M[r = rn^tg/Pg and the mass scale is then 

just m = po,diskr\ = M[r = r//]fg. Unless explicitly stated, the 
magnetic field strength is given in Heaviside-Lorentz units, 
where the Gaussian unit value is obtained by multiplying the 
Heaviside-Lorentz value by \/47r. 

The value of pojisk can be determined for different sys- 
tems. For example, a collapsar model with M = 0. IM05"' and 
M K. 3Mo, then po.d/it w 3.4 x 10'"gcm"^. M87 has a mass 
accretion rate of M ~ I O~^Mp yr~' and a bla ck hole mass 
of M w 3 X lO'^M© (lH^IT99l iR^vnolds et al.llT99^ giving 
p{.),disk ^ 10"'^gcm"^. GRS 1915H-I05 has a mass accretion 
rate of M ^ 7 x IO" ' ^MfT, v r"' (Mirabel & Rodriguez 1991 
'Mirab el & RodrfguezllT999 . Ipender & Belloni 2004) with a 
mass of M ^ 14Mr7^ dGrein er et al. 2001 ), but see Kaise r et alJ 
(|2004j). This gives po.i/Kt 3 x lO'^^gcm"^^. This disk density 
scales many of the results of the paper. 

3. BREAKDOWN OF IDEAL-MHD 

Pair creation is critical to understand the physics of the 
highly magnetized, evacuated funnel region that is associ- 
ated with a Poynting-dominated jet. Pair creation is often in- 
voked in order to use the force-free electrodynamics or ideal 
MHD approximation in a black hole magnetosphere (see, e.g., 
iBlandford & Zna iek 1977). However, in MHD where rest- 
mass is treated explicitly, pair creation is not simply a passive 
mechanism to short out spark gaps, which is the mechanism 
invoked to allow the use of the force-free approximation. 

Pair creation (and neutron Fick diffusion for GRB-type 
systems) determines the matter flow in the magnetosphere, 
and th us the matter-loading o f any Poynting jet that emerges 
( Phin nevI 1 983t iPunslvll 1 99 U ILevinson 2005). As shown be- 
low, these sources of mass loading self-consistently determine 
the Lorentz factor of the Poynting-dominated jet and allows 
one to understand why black hole accretion systems, while 
following the mass-invariant GRMHD equations of motion, 
show a variety of jet Lorentz factors. 

For GRBs, the radiative physics and neutron diffusion 
is shown to determine the Lorentz factor of the Poynting- 
dominated jet. For AGN and x-ray binaries, the radiative 
physics is shown to determine the Lorentz factor of the 
Poynting-dominated jet by determining its energy and mass- 
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loading. 

The ideal MHD approximation (or force-free approxima- 
tion in magnetically dominated regions) has been shown to 
be a reasonably valid theoretical framework to describe most 
of the nonradiative dynamically important accretion physics 
around a black hole i n GRBs , AGN, and black hole x-ray bi- 
nary systems (Phinn e\ll983t [Mc Kinney 2004). This approx- 
imation is the foundation of most studies of jets and winds. 
The ideal MHD approximation is a good approximation to de- 
scribe these flow properties except 1) in current sheets, which 
is not treated explicitly in this paper ; 2) where pair creation 
contributes a nonnegligible amount of rest-mass, internal en- 
ergy, or momentum density; 3) if the Goldreich-Julian (GJ) 
charge density is larger than the number density of charge car- 
riers ; and 4) if the rest-mass flux due to ambipolar and Fick 
diffusion is negligible. 

The first goal is to show that radiative annihilation into pairs 
establishes a density of pairs much larger than the Goldreich- 
Julian density. The Goldreich-Julian charge density is never 
reached because pair creation is completely dominated by 
neutrino annihilation in GRB-type systems and photon anni- 
hilation in AGN and x-ray binary systems. 

Notice that the breakdown of ideal MHD is required in or- 
der to extract black hole spin e nergy from a stationary, ax- 
isymmetric system. 191 A) showed that a rotating black 

hole induces a parallel electric current in the surrounding 
magnetosphere such that the plasma bec omes nondegenerate 
(i.e. E'BiiQ). Bekenstein & OronI ( fT97 8) argued that if the 
ideal MHD approximation were valid, that no energy could 
be extracted from a black hole. This is because since m'' < 
at the horizon, and the radial energy flux can be written as 
-T[ = Epu'' (where E is conserved along each flow line ; see 
appendix IH}, then to extract net energy (-r/ > 0) from the 
black hole requires E <Q. However, in the ideal MHD ap- 
proximation £■ > at r ^ oo, and by conservation of E along 
each flow line, then £ > on the horizon as well. How- 
ever, based upon ar guments by Goldr eich & Jul ian ( 1969), 
iBlandford & Znaiekl(tl97 7) argued that as the magnetosphere 
is evacuated to the Goldreich-Julian charge density, the par- 
allel electric current separates the charges. The Goldreich- 
Julian rest-mass density for a species of electrons is 



PCJ 



Incq ' 



(1) 



where B is the magnetic field strength and q is the electron 
charge. Once the parallel electric current is sufficiently large, 
electrons are accelerated across the potential gap and pho- 
tons can be emitted by curvature radiation or inverse Compton 
scattering. These high energy photons either self-interact or 
are involved in a magnetic bremsstrahlung interaction, ulti- 
mately leading to electron-positron pairs. These pairs would 
continuously short the induced potential difference. However, 
this picture does not establish how the resulting pair plasma 
flow behaves. 

Why have ideal GRMHD numerical mod- 
els demonstrated the Blandford-Z najek ef- 
fect ( Koide, Shibata. Kudoh. & Mei ej 12002 ; 
iMcKinnev & Gammie 20Q4; De VilUers et al.1 12005a : 
iKomissarov 2005) ? These ideal GRMHD numerical models 
implicitly break the ideal MHD approximation in the required 
way to allow the extraction of energy from the black hole. 
Fo r all the initial co nditions and field geometries explored 
bv IMcKinnev & Gammie ( 2004) using the "ideal" GRMHD 
numerical model of an accreting black hole, they always 



find that a highly magnetized polar region forms and any 
material in this magnetosphere is either rapidly driven into 
the black hole or driven out in a wind or jet. They find that 
strong field lines tie the black hole horizon to large radii. 
Thus, necessarily these ideal MHD models break the ideal 
MHD approximation at a stagnation point where the poloidal 
velocity u'' = 0. Necessarily matter is created (at least) in this 
location since matter inside this surface goes into the black 
hole and matter beyond it goes away from the black hole. This 
aspect is similar to the charge-starved magnetosphere models 
where there is a spark-gap (Ruderman & Sutherl and 197^ 
where particles are generated (for a review see iLevinsonI 
I2005h . Once the magnetosphere reaches an axisymmetric, 
quasi-stationary state, then the departure from the ideal-MHD 
condition can be measured as deviations from conservation 
of the conserved flow quantities given in equations IB3I to 

equations 1^21 

Notice that for a realistic accretion disk the BZ power is 
different than the typically used estimates (McKinnev 2005a). 
For j > 0.5, they find that the efficiency in terms of the mass 
accretion rate is 



and 
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15% 



;7% 



^H[j=l] 



(2) 



(3) 



where r^ = GM/c^, flu = jc/ilru) is the rotation frequency 
of the hole, rn = rg{l + j^) is the radius of the horizon 
for angular momentum J = jGM^/c, and j = a/M is the di- 
mensionless Kerr parameter, where -1 < j < 1. However, net 
electromagnetic energy is not extracted for j < 0.5 (includ- 
ing ret rograde a ccretion) when an accretion disk is present 
(McKinnev & Gammie 2004) . This high efficiency is a re- 
sult of the near equipartition of the magnetic field strength 
((B'')^) in the polar region at the horizon and the rest-mass 
density in the disk at the horizon. If the black hole has j w 0.9, 
then w 1 % of the accreted rest-mass energy is emitted back as 
Poynting flux in the form of a jet and w 3% is emitted back in 
total (so obviously 2% goes into the disk and corona - about 
equally it turns out). 

3.1. ORB Pair Creation Model 

In GRB models, such as the collapsar model, neutrino/anti- 
neutrino annihilation provides a source of electron-positron 
pairs at a much larger density than the Goldreich-Julian den- 
sity and so the magnetosphere is not charge starved. The 
cross-field magnetic diffusion for charged species is negligi- 
ble in such systems. However, free neutrons diffuse across 
the field lines and load the jet with an electron-proton plasma 
dLevmson & Eichleil200l . and this effect is considered in the 
next section. 

For the collaps ar model, the jet has B ^ 3 x lO'^Gauss 
( lMcKinnevl2004h and j ^ 0.9, which gives 



Pcj ■ 



lQ-''gcm-^ 



(CoUapsars). 



(4) 



One can compare this to the density of pairs produced by 
neutrino annihilation for the GRB collapsar model. One can 
use the results in table 3 and figure 9b of Popham et al. ( 199^ 
and the results in table 1 in .MacFadven & Woosley ( 1999t) , 
which are fairly well fit to power laws, such that for models 
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with M < 0. 1Mq5"' , and M = 3Mq 



'1% 



M 



(5) 



where this fit is based on an average between the conserva- 
tive and optimistic models of MacFadven & Wooslev ( 1999). 
This assumes an average neutrino energy of ^ lOMeV from 
the disk. This says that for the collapsar model with j = 0.9 
that about 1% of the rest-mass accreted is given back as 
positron-electron pairs due to neutrino annihilation, which is 
similar to the Poynting flux from the black hole that goes into 
the jet region as from equation|3] 

Published results of ne utrino annih ilation rates as a func- 
tion of position (Popham et al.lll999l) can be used to obtain 
a preliminary model of pair creation and incorporated into a 
GRMHD model. The details of this preliminary model end up 
not affecting the results, and a more self-consistent model is 
left for future work. The results primarily depend on the over- 
all annihilation luminosity and the basic radial dependence of 
the energy injected as pairs^ 

Figure 9 and table 3 in IPopham et aH (Il999h and table 1 
of MacFadven & Wooslev (1999) can be used to obtain ap- 
proximate radial and height dependent fits of the energy den- 
sity rate of depositing pairs into the jet region. Their figure 9 
shows that the height and radial dependence of the pair annihi- 
lation luminosity per unit distance. These follow approximate 
power laws or exponential laws for j > 0.2. A reasonable fit 
is that 



and 



P[R] 



Q[z] 



(6) 



(7) 



for the luminosity per unit distance. The coefficient is de- 
termined by the total annihilation luminosity {L^j^ann)- As 
in Popham et al.l (il999l) . photon null geodesic transport in 
curved spacetime is neglected such that 



'27rA 



oo /• arctan \H jR\ 



P\R\Q\z\p- ^mBdMr. (8) 
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Figure 6 of iPopham et al.l (Il999h can be used to obtain the 
disk thickness to radius ratio 



2r„ 



2/3 



(9) 



for M = O.IM0S"'. This allows one to determine that A « 
i'!yp,fl«n/191. Thus the energy generation rate can be written 
as 



■P[R]Q[z] (Collapsars), 



(10) 



where Na ~ 191 and we have defined po.disk = M/r^^. How- 
ever, the above H/R assumes the jet fills around the disk. 
Rat her, there is like ly a thick corona between the disk and 
jet iMcKinnev & Gam mie 2004). Motivated by those simula- 
tions and the simulations discussed in this paper, the jet region 
is presumed to exist within 



0:^1.0 



3rg 



-1/3 



(11) 



for r < lOOrg. In that case Na ~ 70. Thus, a significant frac- 
tion of the pairs are absorbed by the corona. However, the 
coro na could also contribute significantly to neutrino produc- 
tion jRamir ez-Ruiz & Socrates 2005). Thus, 70 < A^a < 191 
are reasonable limits. 

Some fraction of the total energy deposited goes into pair 
rest-mass, pair internal energy, radiation, and pair momentum. 
Here fp denotes the fraction turned into lab-frame mass, //, the 
fraction turned into lab-frame internal energy and radiation, 
and /„, the fraction turned into momentum energy. Thus 1 = 
fp+fh+fm- The pair rest-mass density creation rate is defined 
as 



Pc'-e+ _ r 

~JP A 



POJisk 



PO,disk 



(12) 



where P(.-g+ = po ^-^+u' . 

One can obtain a rough density measure in the injection 
region by assuming the characteristic time scale for moving 
the pairs once they have formed is the light crossing time 
at the stagnation surface tstag ~ tg(rstag/ fg). Then po.e-e+ ^ 

PO,e-e+tstag- With po,disk = POJisktstag, 



Pe-e+ 
PQ.disk 



'fp 



POJiskC' 



^ I . (13) 



Notice that many pairs fall into the black hole, so do not 
contribute to the jet rest-mass or energy. Only those pairs 
beyond the stagnation surface survive the gravity of the black 
hole, such that the total annihilation luminosity into the jet is 



Lesc = 27r 



(14) 



which is a similar integral as performed before. However, no- 
tice that particles injected with r < rs,ag, by definition, fall 
into the black hole since they are inside the stagnation sur- 
face where < 0. This is unlike the BZ-power, which in 
steady state is well-defined and conserved through the stag- 
nation surface (Levinson 2005). For r^nj = r^, all the injected 
mass reaches infinity by definition. For rstag = lOrg, only 11% 
of the mass injected reach infinity. Because any mass injected 
lost to the black hole is of no consequence to the acceleration 
at large distances, then the true efficiency of pairs that load 
the jet is 

^-^ = M? (15) 

rather than rj^i^^am for all rsiag. One can show that rjesc = iliyi^,ann 
for r^flg = rg, but is reduced to r]esc ~ 0.1 l77i,p_n„„ for rstag = 
lOrg due to the loss of pairs into the black hole. 

The pairs annihilate after formed by neutrino annihila- 
tion. Equation ID7I gives the pair annihilation rate. For 
GRB models, such as the collapsar model, the pair annihi- 
lation timescale is tpa ~ 10~'^s ^ GM/c^ ~ 10~^s and sim- 
ilarly all along the jet. Thus, the pairs annihilate and form 
a thermalized fireball. A fraction fp + fh of the energy in- 
jected is turned into a electron-positron pair-radiation fireball. 
The ty pical an gle of scattering neutrinos gives /,„ ^ fp + fh 
(Poph am et alj ri999). Thus, for the fireball formation region 
fp+fh ^ fn ~ 0.5 within factors of a few, and this is indepen- 
dent of the energy of the neutrinos or the efficiency of annihi- 
lation. 

All of the mass energy ther m alize s into the fireball with a 
temperature given by equation IA12I The formation fireball 
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pair mass plus pair internal energy density plus radiation in- 
ternal energy is 

lQ,lotl^'Ut+ Pg-g++ / ^uu.ann \ f ^stag\ (16) 
POMiskC^ \ POMiskC^ J \ ''g J ' 

where qotot = (P(),e-e+c^ + uo,e-e* + uo.^). This equation con- 
nects the energy injection process in terms of the GRMHD 
equations of motion for a given energy-at-infinity injection 
rate. See also appendixlAl 

Using the discussion here and the equations in appendixlAl 
one can show that for a = 0.1, M = O.lM0.r', and j = 0.9, the 
fireball temperature is T ^ lO'^K ^ IMeV in the injection 
region. Thus pair rest-mass energy is nearly in equipartition 
with the pair internal energy and radiation. In particular, fp w 
/,,/8.5. Since fp+f,, - 0.5, then fp « 0.05, « 0.45, and 
f,n - 0.5. 

Beyond the initial fireball formation, the Boltzmann factor 
gives that once the temperature drops below T 6 x lO^K 
the number of pairs decreases exponentially with temperature. 
However, the fireball is optically thin only at much larger radii 
of r ^ 10^- 10"Vg. So until that radius, the radiation provides 
an inertial drag on the remaining pair plasma and the gas is 
radiation dominated. 

For example, for j = 0.9, M = O.IMqs"', rsiag = 6rg, then 
the initial fireball rest-mass is 

-^-10-^ (r,,,g = 6r,). (17) 

Pl).disk 

If instead the model were an a = 0.01 model, then the ef- 
ficiency would be about 10 times less and the density ratio 
would be about 10 times less at 

-^^-10"^ (a = 0.01). (18) 

PO,disk 

In any case this is roughly 

Pe-e- > lO'gCm-\ (19) 

which is about 12 orders of magnitude larger than the GJ den- 
sity given in equation 0] and so the black hole is far from 
starved of charges. 

3.2. GRB Baryon Contamination 

Neutron diffusion across field lines leads to baryon contam- 
ination of the (otherwise) electron-positron-radiation jet. The 
neutrons Fick-diffuse (Levinson & Eichler 2003) or diffusion 
due to ambipolar diffusion (see appendix |XJ across the field 
lines and baryon-load the jet. The neutrons undergo a fast col- 
lisional avalanche into protons and electrons that are then car- 
ried along with the electromagnetic and Compton-thick out- 
flow. As shown in appendix|A| the mass injection rate of neu- 
trons (and so proton-nelectrons) is 

M in j, Pick ^l-X lO'^Macc (20) 

where the density of electron-proton plasma is 

Pp.- - 3 X 10"^ (^^j pQji^k- (21) 

Notice that this mass injection rate and density are com- 
parable to the injection-region rest-mass density in electron- 
positron pairs for a = 0.01. As mentioned above, father out 
in the jet the electron-positron pairs annihilate and contribute 
only an additional ~ 10% to the internal energy. Thus the total 



internal energy is sufficient to describe the gas without includ- 
ing pair annihilation and the rest-mass in baryons is sufficient 
to describe the gas rest-mass for models with a = 0.01 and for 
all models for r > lOr^ where the pair-density is exponentially 
smaller than the baryon density in the jet. Thus, the injec- 
tion of "pairs" described in the previous section can also ap- 
proximately account for the Fick diffusion of neu trons. This 
fact is exploited to simulate the collapsar model in lMcKinnevI 
(2005b). 

3.3. AGN and X-ray Binaries Pair Creation Model 

In AGN and x-ray binaries, the pairs are produced by scat- 
tering of > IMeV 7-rays with other photons (for a review 
see Phinney 1983 and chpt. 6, 9, and 10 in Punsly 2001). 
These 7-rays could be produced by, for example, Comp- 
tonization of disk photons through a gas of relativistic elec- 
trons, non-thermal par ticle acceleration in shocks (see, e.g., 
'Nishi kawa et al.n 2003'). or reconnection events. For exam- 
ple, the (non-radiative) simulations of McKinnev & Gammi^ 
(2004) show an extended corona that could be source of 
Comptonization. They also find an edge between the corona 
and funnel that contains frequent shocks with sound Mach 
number M ^ 100. Also, they found reconnection is quite vig- 
orous in the plunging region at r ^ 3-6rg, leading in part to 
the hot coronal outflow. 

These sites of Comptonization, shocks, and reconnection 
are likely sources of the requisite 7-rays. In place of a detailed 
model of these processes, it is assumed that some fraction (/y) 
of the true bolometric luminosity is in the form of these 7-rays 
that do collide with softer photons in the funnel region. For 
a bolometric luminosity Lhoi ~ rjef/Mc^, then the annihilation 
efficiency is 

Vll.ann = fl j^^2 ~ fl^^ff^ (22) 

where rjeff is the total radiative efficiency, which could in- 
clude emission from both the disk and jet near the base. No- 
tice that rjfff depends on the black hole spin, among other 
things. However, the value of rj^fj is obtained from (model- 
dependent) values for the mass accretion rate and bolometric 
luminosity. 

Extrapolating from gamma-ray observations of black hole 
x-ray binary systems suggests that in either the quiescent or 
outburst phase, ^ 1% of the true bolom etric luminosity is 
in the form of > IMeV photons (see, e.g., iLing & WheatonI 
I2005h . These are likely produced quite close to the black 
hole. For an injection regio n with a half opening angle 
for the jet of 9j ^ 57° (M cKinnev & Gammie 2004). about 
(29j)^/(4TT) w 1/3 of these photons enter the jet region. Thus, 
it is assumed that a large fraction of these > IMeV photons 
give up their energy into producing pairs in the funnel region 
with some fraction of the energy going into rest-mass (fp). 
These pairs do not annihilate so form a collisionless plasma 
(see appendixlAl. See alsojPunslv (200lJ) for why fy ^ 1% 
is reasonable, based on assuming the infall rate is equal to the 
pair creation rate. In general depends on the state of the 
accretion flow, and a self-consistent determination is left for 
future work. 

As in the collapsar case this gives us a density rate or a 
typical density. In this case the stagnation surface is close to 
the black hole since the emission is likely always optically 
thin, thus rs,ag = 3rg. The author knows of no calculations 
that give a radial dependence for the annihilation energy gen- 
eration rate. A reasonable choice is the same radial depen- 
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dence as for neutrino annihilation, since while the efficiency 
of neutrino scattering is much lower, the radial structure is 
determined by a similar calculation. The typical thick disk 
ADAF model assumes the electrons and protons are weakly 
coupled, which gives a disk thickness mostly independent of 
radius such that H /R ^ 0.85 and is weakly dependent on a 
(jNaravan & Yi 1995). For the AGN and x-ray binaries, jets 
are presumed to occur in the presence of a thick (ADAF-like) 
disk close to the black hole. While observations show that the 
disk at r > 6rg undergoes state transitions, steady jets are only 
observed in the low-hard state when the disk is likely geomet- 
rically thick. For this H /R, equation |8| gives A = L-y-,, a„„/53 
or A'a = 53. Now equations 1121 and [Ts] also apply for photon 
annihilation but with Na = 5 3. __ 

For AGN accretion disks, Phinnev ( 1983) already showed 
that 7 rays from the accretion disk corona interact with x-rays 
to produce electron-positron pairs in sufficient density above 
the Goldreich-Julian density. For example, M87 has a nuclear 
bolometric luminosity of L-y « 2 x IC^^erg s"' and a mass ac- 
cretion rate of M lO'^Mgyr"' ("Ho"T99p, iRevnolds et al.1 
[T996) giving rjeff ~ 4 x 10"^. Unlike neutrinos, 7-rays are ef- 
ficient at creating pairs and most of the 7-rays are at around 
IMeV so the fraction of the energy put into rest-mass is 
/p~l.If/^~l%,then 

p,-,. - 10-V'«"' (/7~1%)- (23) 
The fi eld strength in M87 is B 0.1 to 50Gauss (iMcKinnevI 
120041) . For ; = 0.9 this gives that 

PGy<10-V'w"' (M87). (24) 

This is about 9 orders of magnitude lower than the pair cre- 
ation established density, so the black hole is not charge 
starved. For M87, 

- 10"^ (M87). (25) 

For black hole x-ray binaries a similar calculation is 
performed. For example, GRS 1915H-105 has a mass 
accretion rate of M ~ 7 x 10 ~^Mf7^yr~' and a bolomet- 
ric luminosity L ^ lO'^^erg s"' (Mirabel & Rodriguez 1994' 
iMrabel & Rodriguez 1999; Fender & Belloni 2004) with a 
mass o f M ^ HM© ijGreiner et al. 2001), but see lKaiser et al.l 
J2004h . This gives 77,// « 0.26. If - 1% and /p ~ 1, then 

Pe-e- ^ lQ-\cm-^ (/7~1%)- (26) 

The field strengt h is ~ lO^Gaus s if the disk is in the thick 
ADAF-hke state jMcKinnevl2004h . For ; = 0.9 this gives that 

pcy - 10"'V'«"^ (GRS1915+105). (27) 

This is about 9-10 orders of magnitude larger than the GJ 
density, so the black hole is not charge starved. For GRS 
1915H-105, 

-^^-10"^ (GRS1915+105). (28) 

PQ.disk 

Notice that GRS 1915H-105, and many x-ray binaries, are 
more radiatively efficient than most AGN. This means x-ray 
binaries have jets loaded with more pair-mass density per unit 
disk density. This will impact the presence and speed of any 
Poynting-lepton jet, as described in the next section. 

Finally, clearly temporal variations in the disk structure and 
mass accretion rate directly affect the actual pair creation rate 
in the jet region. A self-consistent treatment of this (time- 
dependent) radiative physics is left for future work. 



4. JET LORENTZ FACTOR 

The Lorentz factor of the jet can be measured either as the 
current time-dependent value, or, using information about the 
GRMHD system of equations, one can estimate the Lorentz 
factor at large radii from fluid quantities at small radii. The 
Lorentz factor as measured by a static observer at infinity is 

T = u' = u'^, (29) 

in Boyer-Lindquist coordinates, where no static observers 
exist inside the ergosphere. This is as opposed io W = 
u' ^J-l /g", which is the Lorentz factor as measured by the 
normal observer as used by most numerical relativists. 

For a GRMHD model, t o dete rm ine th e Lorentz factor at 
r ~ 00, notice that eauations IB 1 1 l and lB 1 2l sDecifv that E and L 
are conserved along each flow line. For flows with a magnetic 
field that is radially asymptotically small compared to near the 
black hole, as r ^ 00. Also, for a nonradiative fluid, 

the internal energy is lost to kinetic energy, and so /1 ^ 1 as 
r ^ 00. Thus Too = -Ui{r = 00]. Now, since $ and Vlp are 
conserved along a flow line, then trivially 

Too=E = -hu, + 'i>D.FB^ (30) 

ut=L = hu^ + <^>B^, (31) 

where h = (po + Ug + p)/ pa is the specific enthalpy, $ is the 
conserved magnetic flux per unit rest-mass flux, ilf is the 
conserved field rotation frequency, and B^ is the covariant 
toroidal magnetic field. E and L simply represent the con- 
served energy and angular momentum flux per unit rest-mass 
flux. Notice the matter and electromagnetic pieces are sepa- 
rable. 

Since E is the hydrodynamic plus magnetic energy flux per 
unit rest-mass flux, B^ ^ simply represents the conversion 
of Poynting flux into rest-mass flux and h-^ I represents con- 
version of thermal energy into rest-mass flux. These are what 
accelerate the flow. A lternatively stated for the magnetic term, 
from equation IC 141 the fluid is accelerated by the magnetic 
toroidal gradients associated with the Poynting outflow. 

A rough estimate of the Lorentz factor at infinity Too can be 
estimated by assuming all the enthalpy flux or all the Poynt- 
ing flux is lost to rest-mass flux that reach infinity. Then from 
equation lB7l one can break up the matter and electromagnetic 
numerator and average the numerator and denominator sepa- 
rately to obtain 

= rr'+rr' - (^esc^mu) (^) , (32) 

where here "esc" refers to those pairs that escape the black 
hole gravity. In steady-state, only those pairs beyond the stag- 
nation surface escape to feed the jet. 

Before estimating the value of the Lorentz factor for vari- 
ous systems, the toroidal field is shown to be the source of the 
acceleration in ideal MHD. This also allows a probe of the 
structure of the jet, which is not possible in the ab ove a ver- 
ages. In ideal-MHD, E and L in equations IB 1 1 1 and IB 1 21 give 
the Lorentz factor at infinity and the angular momentum per 
particle at infinity. 

The enthalpy is just 

/z = (p + « + p)/p= 1 + 1^^ = 1 + 1^1^ (33) 
3p 3/p 

for a relativistic gas of electron-positron pairs that have fp of 
energy into rest-mass and the rest into thermal or momentum 
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energy. Now the magnetic term is 



pou'(v'^/nF-iy 



(34) 



Notice that E diverges for v"^ = ilf (or u'' = = 0), where 
the ideal MHD approximation breaks down. Thus, Too is also 
determined by the non-ideal MHD physics of particle creation 
in that region. 

For an extended particle creation region, Too depends on 
the mass, enthalpy, and momentu m injected as a f unction 
of radius and in the jet region (Levinson 2005; Punsl^ 
pool). For a narrow {6r <C rg), yet distributed, particle cre- 
ation region, then plausibly ideal MHD is completely reestab- 
lished at a slightly larger radius where v'^ ^ Q,p k, 57/^/2 
(|McKinnev & Gammie 2004). In this case $ w -pflfu' /B'l' 
(or u''/B'' Ri -ilfu'/B'^). Thus, the magnetic piece is 



<^nFB4, : 



B^ 

P 



(35) 



where p = pou' is the lab-frame density. Written in Boyer- 
Lindquist coordinates in an orthonormal basis, then the mag- 
netic piece is 

r(£M)^(^^ (36) 



where fi"^: 



pB"^. Thus, in Boyer-Lindquist coordinates 



1 + 



4(1- 



3/p 



■/p)^(fi^ 



pc^ 



(37) 



for r > fh- Thus the fluid energy at infinity is due to con- 
version of thermal energy and toroidal field energy into ki- 
netic energy. The actual value of Too depends on how nar- 
row is the injection region and the location of the stagnation 
surface. The relevant magnetic field is the toroidal magnetic 
field strength beyond the injection region where ideal MHD 
is mostly reestablished. Notice that a = b^/ po is often used to 
parameterize the Lorentz factor for a magnetically dominated 
flow, while perhaps equation|^is more useful. 

The polar field on the black hole horizon is approximately 
monopolar even for rapidly rotating holes with j < 0.95 
(|McKinnev & Gammie 2004). The monopole field solution 
can then be used, when properly normalized, to give the func- 
tional dependence of B"^ near the black hole. Then one may 
use equation]^ to obtain the approximate terminal Lorentz 
factor First, the black hole emits a Poynting flux given by 
equation 33 in McKinnev & Gammie (2004) with the use of 
their section 2.3.2 for the magnetic field. In an asymptotic ex- 
pansion, which is good to factors of 2 for any j even at r = 3rj, 
then 



B'^ 



-c- 



smt 



Modifications at higher black hole spin are within factors of 
1.5 in magnitude and there is a factor of < 1.5 enhance- 
ment near the polar axis compared to lower spin. In order 
to use this simple BZ monopole, the jet region's "C" co- 
efficient can be jjbtained from GRMHD numerical models 
jMcKinnev & Gamm ie 2004). For a - 0.2-0.4 disk 
model, they find that the normalization constant "C" in the 
BZ monopole solution is 



C w O.ly/pojiskC 



(39) 



(see equations 47-49 in iMcKinnev & Gammi3 12004 . This 
clearly states that the toroidal field in the polar region is nearly 
in equipartition with the rest-mass density in the disk. One 
can plug this into equation]^ but this equation would have 
limited usefulness since the injection region is broad and the 
density is difficult to estimate. However, notice that Too has 
angular structure if the magnetic term dominates since then 
Too c>c sm^9. This is important to jet structure described in 
McKinnev (2005b). Otherwise the result has the same quali- 
tative features as eauationl46l 

4.1. Lorentz Factors in Collapsar Systems 

This section shows that, without invoking super-efficient 
neutrino emission mechanisms, only the Blandford-Znajek 
driven process can drive the flow to the necessary minimal 
Lorentz factor to avoid the compactness problem and un- 
equi vocally generate a GR B. Some GRBs require up to F ^ 
500 ("Li th wick & Sarfl2001 ), so any invoked mechanism must 
be able to explain this. Some observation/models suggest 
some bursts have even F ^ 1000 (ISoderberg & Ramirez-Ruia 
120031) . 

Equation|32]can be written as 



(Vesc + 1lEMJ. 



(40) 



Equation[^gives rjesc, which accounts for pair capture by the 
black hole. Here tjem is given in equation |3l and as gener- 
ally noted before, 77^5^ = 'nvu,ann for rstag = r^, but is reduced 



to rjesc ~ 0.11?7i/p,a„„ for rstag = lOr^ due to the loss of pairs 
into the black hole. Based upon GRMHD numerical models 



studied in lMcKinnevI ( l2005bl) . a likely value is rstag = 5rg, for 
which rjesc « 0.577i,p,a„„. 

GRB-type systems are different than AGN and x-ray binary 
systems, because neutron diffusion baryon-loads the jet. The 
indirect injection of protons and electrons is shown to dom- 
inate the rest-mass in the jet because the electron-positron 
pairs annihilate to negligible rest-mass . Th e Fick diffusion 
mass injection rate is given by eauation iA 1 01 and is 

Minj^Fick-^Ty-lO-^Macc (41) 

(•Levinson & Eichlei'2003), and so M^^c = Mi„j^fick- 

Plugging in the efficiencies for neutrino annihilation (equa- 
tion |5} and Poynting-dominated jet efficiency (equation |3} 
into equation 1401 one finds for M < O.IMqs'^ for ; > 0.1 
that 



140 



+ gEM 



.jet) 



(42) 



where 



5 fu.ann.esc - 



a 



J_ 

0.9 



M 



3.8 



(43) 

whe re for M > 0. IM^:^^ ' , the power of 3.8 sharply levels out 



wiici c iui ivi u.iivii^A , iiic puw 

to (iDi Matteo et al.l2002l) . Also, 



gEM.jel 



1 



(44) 



These ^'s are just the normalized the efficiencies. It turns out 
that for i w 0.9 that the efficiencies are similar for the shown 
normalization of parameters. Thus, one might expect that they 
contribute equally to the energy of the jet. 

If fsiag = fg and a = 0.1, then for the collapsar model with 
M = 0. IM0S"' and j = 0.9, then neutrino annihilation and BZ 
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power are equal and typically give T ^ 100- 140, which is 
sufficient to avoid the compactness problem for typical bursts. 

However, r^a^ « 5rg is more likely. The neutrino-driven 
term then gives Too ~ 70 and is dominated by the BZ-driven 
term that stays at Too ~ 100. This is a result of the loss of pairs 
into the bla ck hole. Notice that the results of ( Popham et al.. 
[1999; Di Matteo et alJl2002h and others do not include Ken- 
geometry to trace null geodesies. While the efficiency for an- 
nihilation should increase due to gravitational focusing, more 
pairs are also absorbed by the black hole. Thus, it is unlikely 
that gravitational focusing helps the neutrino-driven mecha- 
nism. 

Also, a = 0.1 probably considerably overestimates the vis- 
cous dissipation rate of a true MHD flow. Despite typically 
a = 0.1 being used by many authors studying viscous models 
of disks. Stone & Pringle ( 2001 ) showed that a = 0.01 is more 
representative of MHD disks near the black hole. This implies 
that the energy generation rate in the disk and the generation 
rate of neutrinos is lower by about an order of magnitude. For 
normal neutrinos this gives Too ~ 14, insufficient to explain 
GRBs. 

Another problem for the neutrino annihilation mechanism 
is that, with the inclusion of optically thick neutrino trans- 
port, the efficiency of neutrino emission is another few times 
lower for the collapsar model (Di Matteo et al."2002). This 
gives Too ^ 5, clearly a serious problem for the neutrino an- 
nihilation model of GRBs. 

The annihilation efficiency approximately scales with the 
average neutrino energy (Pocham et al. 1999). One must in- 
voke super-efficient neutrinos with an average neutrino en- 
ergy of ^ 210MeV (Ramirez-Ruiz & Socrates 2005) in order 
to obtain a neutrino annihilation power comparable to the BZ 
power However, this is near the peak neutrino energy esti- 
mated to come from a hot corona ( Ramirez-Ruiz & Socrates 
120051) . and the corona is not expected to be the dominant 
source of neutrinos, so the average neutrino energy should be 
smaller 

Notice that choosing j = I only increases the neu- 
trino efficiency, and so Too, by a factor of 2. Such 
a black hole spin is only achievable when H/R < 0.01 
fOammie. Shaciro. & McKinnev 2004), which is not repre- 
sentative of any GRB model (Kohri et al. 2005). 

So, without evoking super-efficient neutrino emission 
mechanisms and unreasonably large neutrino annihilation ef- 
ficiencies based upon only optically thin emission, one must 
turn to the Poynting flux to drive the jet. For j = 0.9 one 
has Too ^ 100, sufficient to avoid the compactness problem 
is many GRBs and only invokes an obtainable black hole 
spin. Unlike the neutrino annihilation case, the BZ efficiency 
has been computed self-consistently from GRMHD numeri- 
cal models of GRB-type disks (McKinney 2005b). 

Just as the neutrino-case has a proposed super-efficient 
mechanism, there are super-efficient magnetic models that 
would increase the terminal Lorentz factor Eventually af- 
ter accreting a large amount of magnetic flux, the magnetic 
pressure dominates th e ram pressure of the accretion flow and 
suspends the flow (Igumenshchev et al. 2003: Narava n et all 
E003). GRMHD numerical models are suggested to not have 
simulated for long enough to see this effect. In this case a 
larger amount of magnetic flux threads the black hole. In this 
magnetically arrested disk (MAD) model, the field strength is 
comparable to the rest-mass density in this case rather than 
only a fraction of it. Then, the BZ efficiency is 100% for j = 1 
and an estimate of the jet efficiency is about 50%, which is 



^10 times larger than previously. Thus potentially F ^ 1000 
is achievable with j = 0.9. Only super-efficient neutrinos with 
average energy 2100MeV can obtain such a jet power. 

At the moment there is an insufficient study of the neu- 
trino annihilation rates for disks that have large optically thick 
regions, so Fqo is not directly estimated for BH-NS or NS- 
NS collisions. However, at higher mass accretion rates the 
neutrino emission efficiency levels off rather than increasing 
(Di Matteo et al. 2002), suggesting these systems should also 
be dominated by Poynting-flux energy like for collapsars. 

In summary, the neutrino-driven mechanism is probably 
dominated by the BZ power. Corrections due to optically 
thick neutrino-transport, a realistic choice for a based upon 
MHD models, and the loss of pairs into the black hole are the 
primary reducing factors compared to previous expectations. 
However, the BZ-effect generates sufficiently energetic GRBs 
while only invoking an obtainable black hole spin. 

Also, while a hydrodynamic jet mixes and can destroy jet 
structure, an electromagnetic jet can internally evolve and 
could have a distribution of Lorentz factors in the jet. This 
is what is found in McKinnev (2005b). Thus, while the elec- 
tromagnetic jet may on average have F ~ 100, the core of 
the jet may have F ~ 1000. Thus, without invoking super- 
efficient mechanisms, only an electromagnetic BZ-driven jet 
can explain all observed/inferred GRB Lorentz factors. 

Based upon the density scaling from the simulations in 
McKinnev (2005b) that are summarized in section and 
based upon appendix |DJ the fireball is optically thick to 
Compton s cattering until r ^ 10^- 10"^rg. Simulations dis- 
cussed in McKinnevI ('2005b') show that Poynting flux is con- 
tinuously converted to heat in shocks and so all the energy 
flux leads to acceleration. Thermal acceleration occurs un- 
til the fireball is optically thin. This acceleration reaches 
F ~ 100- 1000 before the internal shocks generate the non- 
thermal emission. 

4.2. Lorentz Factors in AGN and X-ray Binary Systems 

This section shows that the disk+jet radiative physics is cru- 
cial to determine the Lorenz factor of jets. It is shown that ra- 
diatively inefficient AGN, such as M87, should have jets with 
2 < F < 10, while radiatively efficient systems, such as GRS 
1915H-105, may have jets with F < 2, but they may be Comp- 
ton dragged to nonrelativistic velocities. 

For AGN and x-ray binaries equation|32]can be written as 

1 f Mc^\ 

Too = 1- (riesc + Vem) N— K (45) 

fp \ ^esc / 

for a fraction M^sc = fpLesc of mass that escapes the grav- 
ity of the hole. Equation [21 shows that the mass-energy 
loading of the jet is reduced, compared to the total injected 
mass-energy, due to loss of pairs into the black hole. Since 

rjesc = Lesc/ {Mc^), then 

r^^lfi + !te£) (46) 

Jp \ Vesc J 

Here r\EM is given in equation [3] As noted before, r\esc = 
r]jf,am for rsrcig = rg, but is reduced to -qesc ~ 0.1 l77T-y,a„„ for 
fsrag = lOr^ due to the loss of pairs into the black hole. 

For these systems rs,ag w 3rg, since the disk is likely opti- 
cally thin and emits harder radiation closer to the black hole. 
This gives ?yyM»-s/f?nmj ^ 0.83 for the ADAF model of H/R. 
Here the fraction of energy in rest-mass is fp ^ 1 due to the 
efficiency of photon-photon annihilation for > IMeV photons 



Black Hole Jet Formation: Theory 



11 



off the plentiful softer photons. As discussed before, obser- 
vations suggest that the fraction of bolometric luminosity in 
7-rays is /y ^ 0.01. Since the radiative efficiency should 
also depend on black hole spin, it is assumed a typical sys- 
tem has / ^ 0.9 since that is a plausible equilibrium spin 
iGammie. Shapiro. & McKinnev 2 004). A detailed model of 
the radiative efficiency as a function of black hole spin is left 
for future work. 
Equation|46lthen gives that 

roo^l + (r?.//)"'. (47) 

For M87, r?,// = 0.004 so Toe- = 250, while for GRS 19 15+105 
r]eff = 0.26 so Too =5. 

As cited in the introduction, the observed apparent Lorentz 
factor of AGN jets is Too S 30, while inferred Lorentz factor 
of some AGN jets is Toe S 200, in basic agreement with the 
above estimates. However, in particular for M87 this is rather 
large. Also, the GRS1915+105 estimate is a bit lai-ge. 

The key difference between the collapsar event and AGN or 
x-ray binaries is that in the collapsar case the photon luminos- 
ity (including Compton upscattered by F^) is negligible com- 
pared to t he jet luminosity, but see Ghisellini et al. ilOOOV, 
iLazzati eTa l. (2004). In M87 the bolometric luminosity Lboi 
is almost that of the jet luminosity Pje,. In GRS 1915+105 the 
bolometric luminosity is greater than the jet luminosity. 

The below discussion is a preliminary check on how radia- 
tive processes affect the above results. The results of simu- 
lations from McKinnev (2005b) are invoked in order to ob- 
tain the density and magnetic structure of the Poynting-lepton 
jet as summarized in section l6T2l The simulations also show 
that at r 10^- 10'* r^, any remaining Poynting flux is shock- 
converted into enthalpy flux until they are in equipartition. A 
self-consistent simulation with synchrotron emission would 
then show the continuous loss of Poynting flux until the syn- 
chrotron cooling timescale is longer than the jet propagation 
timescale. This still suggests the jet magnetic field is finally in 
equipartition, but that much of the energy is lost and so cannot 
acceleration the jet. 

The jet, rather than just disk, radiative physics is necessary 
in order to explain why Too = 250 is not achieved in M87 and 
Too = 5 is not achieved in GRS 1915+105. For M87, much 
of the Poynting energy that leads to T^o ~ 250 is converted 
to heat by shocks, which is then lost to nonthermal and some 
thermal synchrotron emission. Thus the Lorentz factor of the 
Poynting-lepton jet achieved by r 10^- lO^ry is likely the 
maximum obtainable. Numerical models of M87 discussed 
in McKinnev (2005b) show that 2 < T < 10. The Poynting- 
lepton jet in GRS 1915+105 is likely destroyed by Compton 
drag or at best Too 2. Other black hole accretion systems 
with different mass accretion rates and radiative efficiencies 
have to be independently checked. The below discussion of 
the disk+jet radiative physics should be considered prelimi- 
nary since a full radiative transport is necessary to obtain a 
completely self-consistent solution. 

4.2.1. Jet Destruction by Bulk Comptonization 

The previous section showed that jets from AGN and x- 
ray binaries survive loading by pair production from 7-ray 
photons, but the produced jet may not survive Compton drag 
(bulk Comptonization) by the relatively soft photons emitted 
by the disk. 

The simulation-based results a re us ed t o de termine the op- 
tical depth as given in equations ID3I and ID5I to compute the 



perpendicular and parallel optical depths to Comptonization. 
For M87 T|| 6 x 10"'' and tj_ 6 x lO"*" at r 5rg (stagna- 
tion surface where jet starts) and tl ~ 4 x 10"^ at r = 120r^. 
For GRS 1915+105, Tji > 1 1 and t^<A (stagnation surface) 
and tl < 3 at r = \20rg. Thus M87 is not Compton dragged 
by the disk photons, while GRS 1915+105 is likely strongly 
Compton dragged by photons that originate near the base of 
the jet and travel up through the jet or across the jet, or by 
synchrotron self-Compton drag. 

A Compton-dragged jet has a limited Lorentz factor that 
reaches an equilibrium between decelerating and accelerating 
radiative processes. Here is it assumed that most of the disk 
emission is at the base of the jet. Then the relevant scenario 
for GRS 1915+105 is the one where all disk seed photons 
that enter the jet are scattered. An isotropic disk luminosity 
Lhoi shining on a conical jet with half-opening angle 9j dumps 
a luminosity of Lseed ~ SjLt,oi/4- into the jet. The photons 
effectively mass-load the jet and an equilibrium Lorentz factor 
is reached, where 

roo-f^)'''-f^V'' (48) 

\2L,eedJ \ VeffOj J 

for a cold beam of electrons (see, e.g.,|Broderick 2004). The 
thermal Lorentz factor is comparable to the bulk Lorentz 
factor, so thermal corrections are not significant. For GRS 
1915+105 this gives a nonrelativistic velocity (F ^ 1) for the 
jet if most of the emission enters the base. Only if most of the 
emission enters far (r > lO^r^) from the base is up to F ~ 2 
possible. Thus is unlikely, so the Poynting-lepton jet that 
forms in radiatively efficient systems, such as GRS 1915+105, 
are Compton dragged to nonrelativistic velocities. Clearly the 
Lorentz factor is sensitive to the disk thickness, the emission 
from the disk, and the structure of the jet. Thus these esti- 
mates should be treated as preliminary. A self-consistent ra- 
diative transfer calculation is left for future work. 

In summary, a radiatively efficient system loads the jet with 
more pairs from the larger number of 7-rays. This sets the 
maximum possible Lorentz factor to be smaller than for ra- 
diatively inefficient systems. For systems with relatively high 
density jets, such as X-ray binaries, the larger radiative effi- 
ciency also leads to an optically thick jet that can be Compton 
dragged. 

4.2.2. Jet Destruction by Pair Annihilation 

Equation ID7I gives the pair annihilation rate. For AGN, 
such as M87, the pair annihilation timescale is tpa ^ 10"s ^ 
GMjc^ ^ lO'^s and for a jet propagation time tje, ^ r/c, a 
lower limit is tpa/tjet ^10^ all along the jet. Thus, most pairs 
do not annihilate. See appendixElon how this affects the fluid 
approximation. See also Ghisellini et al. ( 1992). 

For X-ray binary GRS1915+105, f^^ ~ 2 x lO-'^s > 
GM/c^ X 10"^ s and tpa/tj^, > 2 all along the jet. Thus, 
some nonnegligible fraction of the pairs annihilate. This also 
contributes to the destruction of the Poynting-lepton jet in X- 
ray binary systems since much of this radiative energy is lost 
at r > 150rg where the jet is optically thin along the jet and 
r > 300rj where the jet is optically thin perpendicular to the 
jet. 

4.2.3. Heat Loss by Synchrotron Emission 

For M87 it was estimated that T^o ^ 250, which is in- 
consistent with observations. However, much of the Poynt- 
ing energy is converted to internal energy in shocks induced 
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by toroidal field instabilities ( lMcKinnevll2005bl) . Thus, the 
synchrotron cooling time might be sufficiently fast to re- 
lease this internal energy that would otherwise accelerate the 
flow throug h thermal a cceleration. For numerical models de- 
scribed in McKinnev (2005b) that correspond to M87, an 
equipartition "magnetic fireball" forms between 10^ rg and 
lO^rg. If this energy could be released, then the shocks would 
again resume and all the Poynting and thermal energy would 
be lost. For an equipartition magnetic fireball half the energy 
is thermal, so the jet internal energy is Ug-i,+ ^ 125/9o,e-f+c^ and 
so the thermal Lorentz factor is ^ 125. The synchrotron 
cooling time in the lab frame is 



^ OTT- 



(49) 



This gives that 



1.4 



tsyn ~ 10^ ( - ) (r < 390r^,M87) (50) 



and 



(r>390r„,M87). 



(51) 



For a typical lab frame jet propagation time of tjet ^ r/c. 



0.4 



and 



!f>!l ^ 10- ( - 

t jet 



(r < 390r„,M87) (52) 



(r>390r„,M87). (53) 



Hence, one would not expect synchrotron cooling to take 
much of the internal energy away. 

However, the "magnetic fireball" forms by shock 
heating and electrons are dramatically accelerated in 
such relativistic coUisionless shocks. The shocks gen- 
erate a power law (nonthermal) distribution of elec- 
trons, where much of the energy is carried by high- 
energy electrons (Beaelman et al. 1984; Blandford & Eichler 
[1987; Achterbera et al. 2001; Fender & Maccarone 20031 
iKeshet & Waxman.2005: .Fender et al...2005.) . Typically the 
distribution is A^(£') oc E~^-^^. For a pair plas ma the maximum 
energy is limited by synchrotron losses ("Achterber y et alJ 
1^01), and the resulting synchrotron emission has photon en- 
ergies of £■ ^ ISTeMeV. This gives 7-ray and up to possibly 
TeV emission beamed along the jet, as in blazars. For exam- 
ple, Mrk 421 shows 15 minute variability, which for a mass 
of 1.9 X 1O^M0 would suggest an emission size on the or- 
der of the horizon size (Punch et al. 1992; Gaidos et al. 1996). 
However, relativistic time effects with rft,,;^ ~' 10 place these 
emissions at r ~ 10^'"?' which coincides with the shock-heated 
transfast region discussed in (McKinnev 2005b). 

Thus, a significant portion of the shock-heated internal en- 
ergy should be emitted by shock accelerated electrons and lost 
through the optically thin jet. In the shocks, inverse Compton 
also contributes to emission of high-energy photons and the 
loss of internal energy. Shock-induced population inversions 
may generate cyclotron masers at shock sites and lead to large 
brightness temperatures (Begelman et al. 2005). 

Thus, it is expected that much of the jet is cold with Thuik ~ 
5-10 left over from pre-shoc k magnetic accelera tion. As de- 
scribed in the simulations of lMcKinnevI (|20053), patches of 



slightly faster or slower bulk F are present by r ^ 10^'"g- In 
M87-based models, these range from 2 < F < 10. 
The synchrotron emission angular frequency is 

3TlqB sin a 



(54) 



where sina ~ 1. This gives a characteristic synchrotron fre- 
quency of 



and 



t/,. - 3 X 10'^ I - 



z/^ ~ 3 X 10'* I - 



-0.7 



Hz (r<390r„,M87) (55) 



-1.5 



Hz (r<390r„,M87). (56) 



For r ^ lO^rg where the fireball begins to form, this gives 
Vc - lOOGHz (radio). By r - lOV^, ~ lOGHz (radio). The 
emission frequency depends on the mass accreti on rate (and 
so po,disk) for any particular AGN. As discussed in lMoKinnevI 
(l2005bl). r ~ IQ^rg is also where the flow goes superfast (su- 
personic). Thus this is consistent with the idea that the radio- 
bright static knots at the base of the jet in, for example, 
Cen A is due to shocks in a transfast (transonic) transition 
(|Hardcastle"'2005). 

In summary, the jet in M87 likely emits most of the internal 
energy, generated in shocks in the transonic transition, as non- 
thermal synchrotron with some thermal synchrotron, such that 
the jet beyond 10^ - 10'*rg is relatively cold with 2 < Too ^10. 

Notice that in x-ray binaries, for example GRS 1915H-105, 
have a jet with u/po^-g+c^ ^ 2.5 and so thermal ~ 2.5. This 
gives that 



^syn 



10-^sf- 



1.4 



(r < 390r„,GRS) 



(57) 



and 



10""s(-) (r>390r^,GRS). 



(58) 



For a typical lab frame jet propagation time of tj,,, ^ r/c, 

OA 



^ 10-^' ( - 

tjet 
tjet 



(r < 390r„,GRS) (59) 



(r>390r„,GRS), (60) 



where GRS denotes GRS1915H-105. Thus thermal syn- 
chrotron is sufficiently fast to cool the jet. Since the jet is 
optically thick, as estimated above, then synchrotron self- 
absorption will dominate the emission process, which is what 
is observed jFosteretalJI19 "96t lFenderfc Rellonill200l . The 
thermal synchrotron emission has 
-0.7 



and 



2x 10' 



j/,.-3x 10" - 



Hz (r <390r^,GRS) (61) 



-1.5 



Hz (r> 390r„,GRS). (62) 



Near the base this gives O.OlkeV emission (EUV). These 
soft synchrotron photons will be Compton upscattered (syn- 
chrotron self-Compton) by the F < 5 jet to x-rays and con- 
tribute to the destruction of the Poynting-lepton jet. Like 
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in AGN, nonthermal synchrotron likely takes away much of 
the shock-generated internal energy and this may account for 
some unidentified EGRET sources. 

It is beyond the scope of the present study to establish 
whether nonthermal synchrotron, synchrotron self-Compton, 
or external Comptonization accounts for most of the high- 
energy luminosity. 

4.2.4. Other issues 

Another possible way of contaminating the Poynting- 
dominated jet is by accreting a complicated field geometry 
and so baryon-loading the polar region. This turns the jet into 
a mixed lepton-baryon Poynting jet. Core-collapse presents 
the black hole with a field geometry that has an overall sin- 
gle poloidal sign. Compared to GRBs, AGN and x-ray binary 
black holes are more likely to accrete nontrivial field geome- 
tries leading to baryon contamination of the jet. Especially in 
Roche-lobe formed disks in x-ray binaries, it is likely that the 
accreted field geometry is quite tangled, so the likelihood of a 
relativistic Poynting jet is further reduced. 

No black hole x-ray binary has been observed to have an 
ultrarelativistic jet (V4641 Sgr is still not confirmed, but see 
[Chatv et al. 2003), despite the GRMHD physics in such sys- 
tems being identical and the Lorentz factor is otherwise inde- 
pendent of the mass of the compact object. However, due to 
their relatively high radiative efficiency compared to AGN, x- 
ray binaries produce more 7-ray flux that increases the pair 
loading for a given magnetic field strength near the black 
hole. Also, the relatively high radiative efficiency means any 
Poynting-leptonjet is severely Compton dragged since the jet 
is optically thick. However, it is possible that there exists 
a large population of low radiative efficiency galactic black 
hole accretion systems. These radiatively inefficient systems 
would produce a large amount of Poynting flux per unit rest- 
mass flux which would be shock-converted by toroidal field 
instabilities into nonthermal synchrotron emission and could 
appear as "microblazars." However, thus-far observed x-ray 
binaries should not be as intrinsically luminous per rest-mass 
accretion rate since the Poynting flux per rest-mass flux avail- 
able to shock-heating is two orders of magnitude smaller than 
available for AGN. This is due to the relatively high pair- 
loading in typical x-ray binaries. Low-luminosity x-ray bi- 
naries would behave more like blazars, and so low luminos- 
ity x-ray binary microblazars may account for some of the 
unidentified EGRET sources. 

5. RELATIVISTIC POYNTING-BARYON JETS 

This section discusses how mildly relativistic Poynting- 
baryon jets can explain many jet observations. The ori- 
gin of these jets is the inner-radial accretion disk. The ori- 
gin of the mass is unstable convective outflows and mag- 
netic buoyancy, and the mass fraction released is typically a 
few p e rcent of the mass accretio n rate (^cKinnev & Gammie 
12004'; 'De Vifliers et al." '2005 a|). The ratio of Poynting to 
baryon flux depends mostly o n the spin of the black hol e 
jPunslv & Coronitil Ll990a b; McKinnev & Gammie 2004- 
Since the Poynting flux from a rapidly rotating black hole 
that is absorbed by the corona is also a few percent, the 
Poynting-baryon jet is heavily baryon-loaded. The heavy 
baryon-loading limits Poynting-baryon jets to only mildly rel- 
ativistic velocities. The most relativistic, collimated, and least 
baryon-loaded portion of the Poynting-baryon jet is at the 
magnetic wall bounded by the Poynting-dominated jet. 



5.1. Matter Jets and Outflows in AGN 

Most AGN should have Poynting-baryon jets. This 
Poynting-baryon jet may often lead to erroneous conclusions 
about the nature of th e jet in AGN systems. 

For example, Juno ret al.l ( fT99 9): Bire ttaet al.l fl999. 200^ 
suggest that M87 slowly coUimates from about 60° near the 
black hole to 10° at large distances. However, two of their 
assumptions are likely too restrictive. First, they assumed 
the jet is always conical, which is apparent from figure 1 in 
Junor et al. ( 1999). If the jet is not conical this can overesti- 
mate the opening angle close to the core (i.e. perhaps 35° is 
reasonable all the way into the core). Second, their beam size 
is relatively large so that factors of 2 error in the collimation 
angle are likely. Finally, and most importantly for this paper, 
they assumed that there is only one jet component. This likely 
leads to a poor interpretation of the observations. If there is a 
highly collimated relativistic Poynting-lepton jet surrounded 
by a weakly collimated Poynting-baryon jet, then this would 
also fit their observations. 

Alternatively, if the accretion dis k in M87 is a ver y thin 
SS-type disk with H/R^ 0.00048 (lMcKinnevll2"()(Ml) . then 
their conclusion that there is slow collimation is plausible. 
However, thin disks may be much less efficient at producing 
jets (Livio, Oailvie, & Prinale 1999; Ghosh & Abramowicj 
1997) and may not be able to produce collimated jets 
( Okam otQ. 1999. .200 ) . A form of the idea that winds coUi- 
mate jets has also been proposed by Tsinganos & Boaovalq^ 
1120051) and applied to M87, but they consider a model where 
the wind slowly collimates the jet in order to fit observa- 
tions. Here we suggest that the observations have been mis- 
interpreted due to the presence of two components: a well- 
collimated relativistic cold Poynting-lepton jet and a mildly 
relativistic coronal outflow. We suggest the broader emission 
component is due to the coronal outflow. 

Notice that more recent maps of the M 87 jet-formation 
region show no "jet formation" structure ( Krichbaum et alJ 
120041) . Thus, the structures seen previously may be transient 
features, such as associated with turbulent accretion disk pro- 
ducing a dynamic coronal outflow. 

Measurements of the apparent jet speed in M87 reveal typ- 
ically r ^ 1.8 near the core while F ~ 6 at larger radii. 
However, so me core regions are associated with F ^ 6 that 
rapidly fade ("Biretta etaDClll)- 

This is consistent with a 
two-component outflow where the cold fast moving core of 
the jet is only observed if it interacts with the surrounding 
medium (or stars), the slower coronal outflow, or it undergoes 
internal shocks. 

For relatively thin disks or slowly rotating black 
holes, Poynting-baryon jets could appear as "aborted jets" 
(Ghiseflini et al. 2004). 

The classical AGN unification models (lUrrv & Padovanil 
fT995h invoke a dominant role for the molecular torus and 
broad-line emitting clouds, while the broad coronal outflow 
may significantly contribute to mo difications and in under- 
standing the origin of the clouds llEivisll20()Ot lElvis et all 

Other erroneous conclusions could be drawn regarding the 
jet composition. Entrainment, which could occur at large 
distances when the ideal MHD approximation breaks down, 
causes difficulties in isolating the "proper" jet component's 
composition. Worse is the fact that there should be two sep- 
arate relativistic jet components, making it di fficult to draw 
clear conclusions regarding the composition (iGuilbert et alJ 
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evinson & BlandfordllT996t 
ISikora & Madeiski 2000). However, is has been recently 
suggested that only electron-positron jets could explain FRII 
sources (Kino & Takahara 2004). 

It is also often assumed that if the jet is highly collimated 
that it is also highly relativistic near the black hole, which 
would suggest Comptonizatio n of disk photons shou ld pro- 
duce clear spectral features tSikora & Madeiskfl2000h . How- 
ever, the jet may rather slowly accelerate and quickly colli- 
mate, which is universally what GRMHD numerical models 
find. 

5.2. Jets in X-ray Binaries 

The results of the previous section suggest that the term 
"microquasar" does not accurately reflect the jet formation 
process. If radiatively efficient systems have no Poynting- 
lepton jet, then what produces their jets? Mildly relativistic 
jets from black hole microquasars may be produced by the 
inner-radial disk rather directly by the black hole. The above 
results suggest that GRS1915H-105 may not have a Poynting- 
lepton jet during its quiescent accretion phase in the low-hard 
state. All black hole accretion systems with a thick disk have a 
mildly relativistic 1 T < 3 coronal outflow due to convective 
instabilities and magnetic buoyancy ( McKinnev & Gammie 
120040 . This component is sufficiently relativistic to explain 
the jets from black hole (and most neutron star) x-ray bi- 
nary systems. This mechanism only requires a thick disk and 
not necessarily a spinning black hole, where other unification 
models suggest that the black hole spin is necessary (Meier 

It has been suggested that the transient, more relativistic, 
jet produced in GRS 1915H-105 is the result the formation of 
a thin disk as t he ADAF collap ses as the mass accretion rate 
increases ( Fen der et al.ll2004bl) . No particular model of the 
transient jet has been suggested. 

Here we give a proposal for the disk-jet coupling in black 
hole x-ray binary systems, such as GRS 1915H-105. In the 
prolonged hard x-ray state the disk is ADAF-like and the 
system produces a Poynting synchrotron self-absorbed jet 
with r < 2, which may be partially or completely Compton 
dragged to nonrelativistic speeds. However, in the thick state, 
a Poynting-baryon jet is prod uced with F ^ 1.5. During the 
soft state, the disk is SS-like ( iShakura & Sunvaevll973 ^ and 
the black hole polar field is relatively weak and the system 
generates an uncollimated (more radial) weak optically thin 
Poynting outflow. During this phase there is also a weak non- 
relativistic uncollimated Poynting-baryon outflow. 

During the transition between hard and soft x-ray states 
the production of pairs decreases significantly in the fun- 
nel, but the black hole polar magnetic field has yet to decay. 
During this transition, an optically thin Poynting-lepton jet 
with r ^ 2 - 3 is produced that is collimated by the remain- 
ing inner-radial ADAF-like structure or the Poynting-baryon 
wind that was produced prior The Lorentz factor produced in 
the transition depends on the details of the disk structure, and 
so 7-ray emission, during the transition. Once the black hole 
field has decayed, the fast transient jet shuts down. 

Alternatively, during the transition to the high-soft state a 
transient jet can emerge as the corona is suddenly exposed to 
more Poynting flux from the black hole. This last bit of coro- 
nal material can be launched off as a faster transient baryon- 
loaded jet. The dynamics of the state transition is left for fu- 
ture work. This overall picture is in basic agreement with 
iFender & Bellonil (|2003), with the additional physics of pair 



creation dominating the Poynting-lepton jet formation pro- 
cess. 

It is interesting that the results ofl Gierliiiski & Don^ ( l2004h 
suggest that for at least some black hole x-ray binaries that 
have jets, the black hole is likely not rapidly rotating (i.e. 
perhaps j < 0.5). For such black holes, there is negligi- 
ble Poynting flux in the form of a Poynting-dominated jet 
dMcKinnev & Gammid 120041) . Thus, our conclusion that 
black hole x-ray binary jets are driven by coron al outflows is 
consistent with the results ofGierli riski & Doli e ( 2004). How- 
ever, even if black hole x-ray binaries were rapidly rotating 
they might not produce Poynting-dominated jets. 

SS443 is plausibly an M 20M(r, black hole system that 
has a jet with v ~ 0.3c (Lopez et al. 2003). Such a low jet ve- 
locity can be explained by a Poynting-baryon jet. To explain 
the opening angle of ~ 1 ° the disk should be very thick near 
the black hole, while the pulsed jet features can be explained 
as an instability due to the overly thick disk self-interacting at 
the poles near the black hole. 

6. SUMMARY OF COMPANION PAPER NUMERICAL RESULTS 

This section summarizes the results of McKinnev (2005^ 
using a GRMHD code HARM (^Gammie et al..200 3a) with an 
advanced inversion method (Noble et alJ l2005ft r Kerr-Schild 
coordinates were used in order to avoid numerical artifacts 
associated with causal interactions between the inner-radial 
boundary and the rest of the flow. Viscous models have found 
this issue to be critical to avoid spurious fluctuations in the 
jet (McKinnev & Gammie 2002), such as might be associated 
with codes using Boyer-Lindquist coordinates. 

6.1. Jet propagation 

As described in detail in lMcKinnevl l l2?)?)5hl) and as shown 
in figure |3l and |3 the Poynting-dominated jet forms as the 
differential rotation of the disk and the frame-dragging of the 
black hole induce a significant toroidal field that launches ma- 
terial away from the black hole by the same force described 
in eauation lC14l 

A coronal outflow is also generated between the disk and 
Poynting-dominated jet. In this model the coronal outflow 
has Foe 1.5. The coronal-funnel boundary contains shocks 
with a sonic Mach number of Ms ^ 100. The inner-radial 
interface between the disk and corona is a site of vigorous 
reconnection due to the magnetic buoyancy and convective 
instabilities present there. These two parts of the corona are 
about 100 times hotter than the bulk of the disk. Thus these 
coronal components are a likely sites for Comptonization and 
nonthermal particle acceleration. 

Figure|3land figurel^show the final log of density and mag- 
netic field projected on the Cartesian z vs. x plane. For the 
purpose s of properly visualizing: the accretion flow and jet, we 
follow MacFadven & Wooslev ( 1999) and show both the neg- 
ative and positive .^-region by duplicating the axisymmetric 
result across the vertical axis. Color represents log(/9o/po,dM) 
with dark red highest and dark blue lowest. The final state 
has a density maximum of po « 2 p^, disk and a minimum of 
/9() ^ po.disk at large radii. Grid zones are not smoothed 
to show grid structure. Outer radial zones are large, but outer 
6 zones are below the resolution of the figure. 

Clearly the jet has pummelled its way through the surround- 
ing medium, which corresponds to the stellar envelope in the 
collapsar model. By the end of the simulation, the field has 
been self-consistently launched in to the funnel region and has 
a regular geometry there. In the disk and at the surface of the 
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Fig. 3. — Jet has pummelled its way through presupernova core and through 
1/1 0th of entire star. Panel (A) shows final distribution of log po on the Carte- 
sian plane. Black hole is located at center Red is highest density and black 
is lowest. Panel (B) shows magnetic field overlayed on top of log of density. 
Outer scale is r = W^GM/c^. 

disk the field is curved on the scale of the disk scale height. 
Within r < lO^r^ the funnel field is ordered and stable due to 
the poloidal field dominance. However, beyond r ^ lO^rg the 
poloidal field is relatively weak compared to the toroidal field 
and the field lines bend and oscillate erratically due to pinch 
instabilities. The radial scale of the oscillations is lO^r^ (but 
up to 10^ rg and as small as lOrg), where r ^ lOr^ is the radius 
where poloidal and toroidal field strengths are equal. By the 
end of the simulation, the jet has only fully evolved to a state 
independent of the initial conditions at r w 5 x lO^r^, beyond 
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Fig. 4. — Strongest magnetic field near black hole in X-configuration due 
to Blandford-Znajek effect and collimation of disk+coronal outflow. As in 
figurel3l but outer scale is r = W^GM/c^. Black hole is black circle at center 
Color scale is same as in figurelsl 



which the jet features are a result of the tail-end of the initial 
launch of the field. The head of the jet has passed beyond the 
outer boundary of r = IQi^rg. Notice that the magnetic field 
near the black hole is in an X-configuration. This is due to the 
BZ-effect having power Pj,,, oc sin^ 6, which vanishes at the 
polar axis. The X-configuration is also related to the fact that 
the disk-i-corona is collimating the Poynting-dominated jet. 
The field is mostly monopolar near the black hole, and such 
field geometries decolUmate for rapidly rotating black holes 
in force-free electrodynamics fcrasnopolskv et alj2005l) . 
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6.2. Summary of Fits 

A summary of the fits along a fiducial field line is given. 
Near the black hole the half-opening angle of the full 
Poynting-dominatedjet is 0j ~ 1.0, while by r ~ 120rg, dj ~ 
0.1. This can be roughly fit by 



-0.4 



(inner) 



(63) 



for r < l2Qrg and 9j ^ 0.14 beyond. The core of the jet fol- 
lows a slightly stronger collimation with 

0j cx (64) 

up to r < 120rg and Oj ^ 0.09 beyond. Also, roughly for M87 
and the collapsar model, the core of the jet has 

^buik ^ ( 57; j (inner) (65) 

for 5 < r < 10^ rg and constant beyond for the M87 model 
if including synchrotron radiation, while the collapsar model 
should continue accelerating and the power law will trun- 
cate when most of the internal and Poynting energy is lost 
to kinetic energy and the jet becomes optically thin at about 
r ~ lO'^r^ or internal shocks take the energy away. If the 
acceleratio n is purely the rmal without any magnetic effect, 
then r cx r jMeszaros & Rees.. l997) . However, it is not clear 
how the equipartition magnetic field affects the acceleration. 
Roughly for GRS 1915H-105 the core of the jet has 

^buik ^ yj^ j (inner) (66) 

for 5 < r ^ ^^^fg constant beyond, with no account for 
Compton drag or pair annihilation. Also, for any jet system 
the base of the jet has po cx r""-' (inner) for r < 120rs, and 
po cx r"^^ (outer) beyond. For the collapsar and M87 models 



and 



PO 



PO 



1.5 X 10" 



1.5 X 10" 



120r, 



-0.9 



-2.2 



(inner) 



(outer), 



(67) 



(68) 



Pojisk Vl20rg^ 

while for GRS 1915H-105 the inner-radial coefficient is 10"^ 
and outer is 6 x 10"^. For the collapsar model, the inner radial 
internal energy density is moderately fit by 

-1.8 



u 



= 4.5 X 10" 



-9 ( 



(inner). (69) 



The outer radial internal energy density is moderately fit by 

-1.3 



u 



:4.5 X 10" 



(outer). (70) 



pi).diskc" V 120rg 

The transition radius is r « 120 rg. For M87 the internal en- 
ergy is near the rest-mass density times until r ^ 120rg 
when the dependence is as for the collapsar case. For 
GRS1915H-105 the internal energy is near the rest-mass den- 
sity times until r ~ 120r^ and then rises to about 2.5 times 
the rest-mass density times c^. The inner radial toroidal lab 
field is well fit by 



V PO,di 



:[Gauss]= 0.0023 



390r, 



-0.7 



(inner) (71) 



for 5 < r < 390rg. The outer radial toroidal lab field is well fit 
by 



V POJiskC- 



: [Gauss] =0.0023 



390r, 



(outer) (72) 



for r > 39Qrg. 

For the typical jet with no atypical pinch instabilities, the 
energy and velocity structure of the jet follow 



e{9) = eoe- 



(73) 



where eo ~ 0.18 and 6*0 ~ 8°. The total luminosity per pole 
is Lj « 0.023Moc^, where 10% of that is in the "core" peak 
Lorentz factor region of the jet within a half-opening angle of 
5°. Also, Too is approximately Gaussian 



^oo(^') = Too.oe 



(74) 



where Too.o ~ 3 x lO-' and w 4.3°. Also, F is approximately 
Gaussian 

r(6i) = roe"^'/^^», (75) 

where Fq « 5 and 0o ~ 11°- The outer sheath's (9 « 0.2) seed 
photon temperature as a function of radius is 



' 50keV 



5 X 10V„ 



-1/3 



(76) 



7. DISCUSSION 



For GRB jets, the injected Poynting flux probably domi- 
nates the injected energy flux provided by neutrino annihi- 
lation. This poses problems for the classic neutrino-driven 
fireball model. Super-efficient neutrino emission mechanisms 
with an average neutrino energy of 210MeV are required in 
order for the neutrino annihilation energy to be as large as the 
energy provided by the BZ effect. However, the BZ effect 
itself might operate in a super-e fficient mode o nce flux has 
accumulated near the black hole JNaravan et aT] l2003'). This 
vertical field threading the disk leads up to 5 times larger lumi- 
nosity (McKinnev & Gammie 2004), in which case an aver- 
age neutrino energy of lOOOMeV is required to compete with 
the BZ-effect. 

For GRBs, equation]^ shows that slightly less rapidly ro- 
tating black holes would generate slightly less Lorentz fac- 
tors but achieve a lower luminosity. This is consistent with 
the observation that harder long-duration bursts have higher 
luminosity, and so suggests that small changes in the stellar 
rotation rate might determine the hardness of long-duration 
bursts. 

The fact that blazars are less luminous for increasing hard- 
ness could be explained by the dependence on the jet Lorentz 
factor on the pair creation physics. Blazars could have similar 
black hole spin, but the hardness of their emission is deter- 
mined by the jet Lorentz factor Lower luminosity systems 
load the jet with less pairs and so the Lorentz factor is larger 
Compton drag of environment or disk reflected seed photons 
can then upscatter to very large energy, such as observed in 
TeV-emitting BL-Lac objects. 

Our results suggest that radiatively efficient x-ray bina- 
ries, such as GRS1915H-105, may only exhibit a relativis- 
tic Poynting-baryon jet. In particular, such a jet is relativis- 
tic only in the low-hard state when the disk is geometrically 
thick. 
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8. CONCLUSIONS 

Primarily two types of relativistic jets form in black hole 
(and perhaps neutron star) systems. The Poynting-dominated 
jet region is composed of field lines that connect the rotating 
black hole to large distances. Since the ideal MHD approxi- 
mation holds very well, the only matter that can cross the field 
lines are neutral particles, such as neutrinos, photons, and free 
neutrons. 

The primary differences between GRBs, AGN, and black 
hole x-ray binaries is the pair-loading of the Poynting- 
dominated jet, a similar mass-loading by free neutrons in 
GRB-type systems, the optical depth of the jet, and the syn- 
chrotron cooling timescale of the jet. 

For GRB-type systems the neutron diffusion flux is suffi- 
ciently large to be dynamically important, but small enough to 
allow r ~ 100- 1000. Beyond r ~ lOr^ many of the electron- 
positron pairs annihilate, so the Poynting-dominated jet is 
dominated in mass by electron-proton pairs from collision- 
induced neutron decay. Most of the energy is provided by the 
BZ effect instead of neutrino-annihilation. 

For AGN and x-ray binaries, the density of electron- 
positron pairs established near the black hole primarily de- 
termines the Lorentz factor at large distances. Radiatively in- 
efficient AGN, such as M87, achieve 2 < Too < 10 and are 
synchrotron cooling limited. The lower the 7-ray radiative 
efficiency of the disk, the more energy per particle is avail- 
able in the shock-zone. Radiatively efficient systems such 
as GRS1915H-105 likely have no Poynting-lepton jet due to 
strong pair-loading and destruction by Comptonization by the 
plentiful soft photons for x-ray binaries with optically thick 
jets. However, all these systems have a mildly relativistic, 
baryon-loaded jet when in the hard-low state when the disk 
is geometrically thick, which can explain jets in most x-ray 
binary systems. 

In an companion paper McKinnev (2005b), a GRMHD 
code, HARM, with pair creation physics was used to evolve 
many black hole accretion disk models. The basic theoretical 
predictions made in this paper that determine the Lorentz fac- 
tor of the jet were numerically confirmed. However, Poynting 
flux is not necessarily directly converted into kinetic energy, 
but rather Poynting flux is first converted into enthalpy flux 
into a "magnetic fireball" due to shock heating. Thus, at large 
distances the acceleration is primarily thermal, but most of 
that thermal energy is provided by shock-conversion of mag- 
netic energy. In GRB systems this magnetic fireball leads to 
thermal acceleration over an extended radial range. The jets 
in AGN and x-ray binaries release this energy as synchrotron 
and inverse Compton emission and so the jet undergoes neg- 
ligible thermal acceleration beyond r ^ 10^- lO^r^. 

Based upon prior numerical (McKinnev 2005b) and this 
theoretical work, basic conclusions for collapsars include: 

1 . Black hole energy, not neutrino energy, typically pow- 
ers GRBs. 

2. Poynting-dominated jets are mostly loaded by e~e^ 
pairs close to the black hole, and by e~p pairs for 

3. BZ-power and neutron diffusion primarily determines 
Lorentz factor 

4. Variability is due to toroidal field instabilities. 

5. Poynting flux is converted into enthalpy flux and leads 
to the formation of a "magnetic fireball." 

6. Patchy jet develops 10^ < Too 10^, as required by 
internal shock model. 



7. Random number of patches (< 1000 for 30 second 
burst) and so random number of pulses. 

8. Energy structure of jet is Gaussian with 6*0 ~ 8°. 

9. Core of jet with 9j w 5° can explain GRBs. 

10. Extended slower jet component with 9j « 25° can ex- 
plain x-ray flashes. 

1 1 . Coronal outflows with F ^ 1 .5 may power supernovae 
(by producing, e.g., ^^Ni) with M ^ O.IM© processed 
by corona. 

Based upon prior numerical (lMcKinnevll20()5bll and this 
theoretical work, basic conclusions for AGN or x-ray binaries 
include: 

1. Poynting-dominated jets pair- loaded unless advect 
complicated field. 

2. 7-ray radiative efficiency, and so pair-loading, deter- 
mines maximum possible Lorentz factor 

3. Poynting-lepton jet is collimated with 6 j w 5°. 

4. Extended slow jet component with 9j < 25°. 

5. For fixed accretion rate, variability is due to toroidal 
field instabilities. 

6. Poynting flux is shock-converted into enthalpy flux. 

7. In some AGN, shock heat in transonic transition lost 
to synchrotron emission and limits achievable Lorentz 
factor to 2 < r < 10 (e.g. in M87). 

8. Coronal outflows produce broad inner-radial jet fea- 
tures in AGN together with well-coUimated jet compo- 
nent (e.g. in M87). 

9. In some x-ray binaries, Compton drag loads Poynting- 
lepton jets and limits Poynting-lepton jet to F < 2 or jet 
destroyed. 

10. In some x-ray binaries, Poynting-lepton jet optically 
thick and emits self-absorbed synchrotron. 

1 1 . Coronal outflows have collimated edge with F < 1 .5. 

12. Coronal outflows may explain all mildly relativistic and 
nonrelativistic jets in radiatively efficient systems (most 
x-ray binaries). 

For AGN and X-ray binaries, the coronal outflow collimation 
angle is strongly determined by the disk thickness. The above 
assumed H/R ^ 0.2 near the black hole and H/R ^ 0.6 far 
from the black hole, while ////? ^ 0.9 (ADAF-like) is perhaps 
more appropriate for some systems. The sensitivity of these 
results to ////? is left for future work. 
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APPENDIX 
A. PAIR CREATION NOTES 

The electron-positron pair plasma that forms may annihilate itself into a fireball if the pair annihilation rate is faster than the 
typical rate of the jet (c^ / GM) near the black hole. Also, if the pair annihilation timescale is shorter than the dynamical time, then 
pair annihilation would give a coUisional term in the Boltzmann equation. From the pair annihilation rate given by eauation lD7l 
one finds that tpa 3> GM/c^ for AGN and marginally so for x-ray binaries. Thus, pairs mostly do not annihilate, and so formally 
the pair plasma that forms in the low-density funnel region is collisionless so that the Boltzmann equation should be solved 
directly. Plasma instabilities and relativistic collisionless shocks are implicitly assumed to keep the pairs in thermal equihbrium 
so the fluid approximation remains mostly valid, as is a good approximation for the solar wind (see, e.g. iFeldman & MarschI 
[T997: Usmanov et al. 2000). This same approximation has to be invoked for the thick disk state in AGN and x-ray binaries, such 
as for the ADAF model (McKinnev 2004). For regions that pair produce slower than the jet dynamical time, each pair-filled 
fluid element has a temperature distribution that gives an equation of state with P = po^e-e+khTe/me rather than P = (ll/12)flr^, 
where a is the radiation constant. So most of the particles have a Lorentz factor of m /(po.e-e+c^) and little of the internal 
energy injected is put into radiation. This also allow the use of a single-component approximation. A self-consistent Boltzmann 
transport solution is left for future work. 

On the contrary for GRB systems, due to the relatively high density of pairs, the time scale for pair annihilation is tpa ^ GM/c^ 
along the entire length of the jet. Thus a pair fireball forms and the appropriate equation of state is that of an electron-positron- 
radiation fireball. Thus, formally the pair fireball rest-mass density is not independent of the pair fireball internal energy density. 
However, because the pairs are well-coupled to the radiation until a much larger radius of r ^ 10*^ - lO'^r^, the radiation provides 
an inertial drag on the remaining pair plasma. That is, the relativistic fluid energy-momentum equation is still accurate. So the 
effective rest-mass density is ^ pu + u (u the total internal energy of the fireball), and so the e ffective rest-mass is independent of 
the cooling of the fireball until the fireball is optically thin (see, e.g.. lMeszaros & Reesll997l) . 

For GRB systems, the mass conservation equation is reasonably accurate. Even though the electron-positron pairs annihilate, 
the rest-mass of pairs injected is approximately that of the pairs that are injected due to Fick-diffusion of neutrons (see next 
section). The annihilation energy from electron-positron pairs contributes a negligible additional amount of internal energy, so 
can be neglected, especially compared to the Poynting energy flux that emerges from the black hole. Thus, the rest-mass can 
always be assumed to be due to baryons rather than the electron-positron pairs. This also suggests that the neutrino annihilation 
is a negligible effect if the BZ power is larger than the neutrino annih ilation power. 

In summary, the rest-mass evolution discussed in McKinnev (2005b) are accurate for GRB, AGN, and marginally so for x-ray 
binaries. This is despite the lack of Boltzmann transport for the collisional system, or a coUisional term due to pair annihilation. 

A.l. Baryon Contamination 

Notice that some frac tion of baryons contamina te the jet due to neutron diffusion and subsequent collisional cascade into an 
electron-proton plasma (Levinson & Eichler 2003). They estimate the diffusion using Fick's law. First, the role of ambipolar 
diffusion is considered (see, e.g. iShu.l992t chpt. 27). The drift velocity is 

Vdrift.pn^-. (Al) 



where L ^ r{H/R) is the typical field radius of curvature induced by disk turbulence and 



nip + p/nin 

is the drag coefficient and (crv) ~ 40c x 10"^^. The drift velocity can also be written as 

Assuming all the diffused neutrons are converted to protonsH-electrons and carried with the outflow, then the diffusion flux is 

111 

F = pVdrif,,pn = — I \" ■ (A4) 

pc^ {(7v)r(H /R) 



and the mass flux across an axisymmetric conical outflow is 



I" out 



Minj.ambi = ^T^ / Frdr (A5) 
Jo 

and so 

Min i,ambi = 2TT{Fr)rout (A6) 

GRB numerical GRMHD models show that the coronal region next to the Poynting- dominated jet has a time-averaged value 
of K pc^) ^ 1 and the turbulent induced eddies occur when the disk has H /R ^ 0.2 (iMcKinnev & Gammiell2004l iMcKinnevI 
12004ft . This gives a mass flux vs. radius of 

MinjMmhi ~ 10"'^ (''_Zjjm\ M^^^^ (A7) 
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where parrs that enter insid e r < rsiag fall into the black hole and do not load the jet. For a recombination radius of r„ ~ 2 x lO^cm 
iLevinson & Eichleil2003h . this gives 

Minj,ambi ~ 4 X IQ-'^Macc (A8) 

This can be compared to the result of iLevinson & Eichled J2003h for Fick-diffusion, where the mass injection rate of free- 
streaming particles (thek eq. 7) is 

Minj.Fick - 3 X 10-' j M„,, (A9) 

for their collapsar model with jet half-opening angle of Oj = 0.1, specific enthalpy h ^ \, Lj ~ 10^' erg s"', a neutron thermal 
velocity of v ~ 0.1c. The thermal velocity is based upon the near-funnel coronal value of u/{poc^) ~ 0.01-0.1 as measured from 
GRMHD numerical models. Notice that this ratio is typically 0.01 in the corona, but is 0.1 at the edge, so we use 0.1 since the 
Fick diffusion is based upon the edge values. Notice that they used v ^ c. For a recombination radius of r„ ^ 2 x lO^cm after 
which no more free neutrons exist, one has that 

7 X lO-^M^,,. (MO) 

Hence, Fick diffusion dominates ambipolar diffusion. The role of reconnection, between the corona/coronal wind and jet, in 
loading the jet with baryons is left for future work. 

The characteristic timescale for moving these pairs is ~ tg{r/rg) and characteristic length is ~ rg(r/rg), so a characteristic 
density vs. radius is 

Mi„i,Ficktg ( ^ ,^-1 ( r\'^'^ 

Ppe ^3 y — j 3 X 10 ' — j pildisk- (All) 

Notice that this is comparable to the rest-mass in pairs given by equation^] Thus, as the fireball decays in pair rest-mass, the 
rest-mass quickly becomes dominated by neutrons diffusing across the magnetic wall between the corona and funnel. Hence, the 
baryon conservation law holds and the approximations used here hold well. A pair-annihilation term is only needed to account 
for the contribution to the interna l energy. Since ^ 8//,, this contribution is a ~ 10% effect and is not expected to affect the 
results of the numerical models of iMcKinnevI i2005bl) . 

A. 2. Pair-Radiation Equation of State 
The total amount of comoving energy put into the thermal fireball is 



AT'' / dx (A12) 
Jo eV'-=+'«- + 1 

where mo,7 = 1.62348Ar'*, A = 4.66244 x lO-'^ergcm'-^K"^, x = pc/khT, m = nieC^/kbT, p is the momentum in the fluid frame, 
and the rest-number density of photons is n^ = 20.2944cm"^r^. The rest-mass in pairs is 

f°° x^ 

po,e-e*=BT^ dx^= , (A13) 

where B = 3.07589 x lO'^^gcm'-'K"-^, and the number density of pairs is no ,,-e+ = PQ.e-e+ /nte- Notice that the GRMHD equations 
of motion relate the comoving energy to energy at infinity by 

Utot = Uo,toiu'u, + Pgas = (fp+ fh + fm)e„i},ann (A14) 

and pe-e+ = pu.e-e+u', where 



AT^ r°° ^ x^V^T^ ^^^^^ 
Pgas=^ dx—==^—+po^^, (A15) 

and p^ = Mo,-y/3. This gives sufficient information to solve for u' and m, or pe-e* = d/dt(pQ,e-e*u'). The study of lMcKinnevH2005bl) 
uses a fixed 7-law gas equation of state with 7 = 4/3 to model the typically radiation-dominated system. 

B. IDEAL MHD QUANTITIES CONSERVED ALONG EACH FLOW LINE 

Kerr spacetime is stationary and axisymmetric with 2 Killing vectors C/^ = ^ = -5'^ and ifi = -§^ = ^'^ that satisfy /^^(g) = 0, 
where C is the Lie derivative and g is the metric. For a vector X^, tensors T that obey £x(T) = are conserved along X. In 
particular, for X = ^, such a tensor is a physical quantity independent of the ignorable coordinates t and 0. For X = u, the 
4-velocity, the tensor is cons erved along each flow line . One can derive a set of conserved flow quantities that are associated 
with the Killing symmetries ( Bekenstein & Oronlll978^ . The below summarizes those results that are key to this paper This 
presentation is necessary for the discussion regarding the determination of the Lorentz factor of the jet. 

In the ideal MHD approximation one can show that £„(C^A^) = 0, where = A^^-A^ defines the vector potential A^. 
Thus for an unsteady axisymmetric flow the (f) component of the magnetic vector potential (A^ = C^^m) conserved along 
each flow line, while for a steady non-axisymmetric flow the electric potential (A, = ^/'A^) is conserved along each flow line. 
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iBekenstein & Oronl (Il978h also show that there are four other independent conserved scalar quantities for an axisymmetric, 
stationary fluid. These correspond to the first integrals of Maxwell's equation f'^^^^ = and the conservation equation T^^ = 
with the use of the continuity equation (/9oM^);p = 0. These first integrals are associated with the projection of the 2 Kilhng vectors 
and the magnetic field {b^^) on these equations of motion. 
For a degenerate ( F F^y = e'^bf^ = 0), stationary and axisymmetric plasma, one finds 

A0,eA,,,-A,.0A0,, = O. (Bl) 

It follows that one may write 

^ = ^^-n,(r,0) (B2) 

wh ere 17^ is usually interp reted as the "rotation frequency" of the electromagnetic field (this is Ferraro's law of isorotation; see 
e.g. lFrank. King. & Rainel2002. §9.7 in a nonrelativistic context). Notice that flp = Ft,-/ F,-^ = Fto / Fee/,- One can show in the ideal 
MHD approximation that 

J7f = y"^ - B^y^ /B^ = y"^ - B'l'v'/B'' (B3) 

is conserved along each flow line. The first term corresponds to fluid rotation and the second term corresponds to the slip along the 
toroidal component of a field line. This yields F^^ in terms of the free functions f7/r,A0, and B'^. Thus, F^ = y/^Er = -ftfA^ ,., 
Fg, = y/^Eg = -il/rA^.e, Fi-e = ^f^B'>', F^r = y/^B^ = -A^j; and Fg^ = y/^B''=A^fi. The diagonal components are zero and 
Fcf,, = £0 = 0, where B' = f' , Ei = Fn/y/^ such that B'Ei = F'^'^Ffj_y/(4y/^^ = E B. Thus for fixed poloidal magnetic field, 
Vlf is a measure of the electric field. With the Faraday written in terms of B' and Vlf, the electromagnetic field automatically 
satisfies the source-free Maxwell equations. 

Using similar constraints on F'"' = A'''*' -A^ '' and with E' = F"^/^, B, = Fj,, one can show that ^J—gF" = E'' = tqBq, 
^pet = -TrBr, ^F't" = E'l' = T^B^, ^f^F''" = -B4,, and ^F'''' = Bg, ^^^^"^ = -Br- Here there is only one 

independent quantity among the three t,'s that are set by F^^F'^'^ = and that the flow and metric are stationary and axisymmetric 
(i.e. E^ = Q). One can show that = BrBgijr - Tg) /B^ and solve for another by using E^ = F°'^gatgi34, = 0. 

It is interesting to note that in Boyer-Lindquist coordinates = and = Tg ^ t and then the contravariant and covariant 
Faraday take on the same simple form with 

where i^zAMO = g''^/g" = 2ar/A is the angular frequency of a zero angular momentum observer (ZAMO) and 

n^^^g^^/g'^ = "-^l^. (B5) 

Notice that if i7f = flzAMO, then r = and so E' = 0. So the difference between the dragging of inertial frames and the field rotation 
frequency generates the electric field E'. Also in Boyer-Lindquist By = B''\g\g^^g't"^(QF/^Q-l), where {glg^^g"^"^ = ('E-2r)/A. 
Also, Bg = B''\g\g''g^H^F/no-l), where \g\g''g'^'^ = S-2r. Also, B^ = -B^\g\g''g'>' , where \g\ = |Det(g^,)| and |g|gV' = 
Asin^e. 

Using the definition of F^'^ given above one can define A^ in terms of the poloidal B'' and B^ giving 

A4lf)-A4li) = {^{B'-d0-B''dr)) (B6) 

over the line segment from /, to / /. Thus given the poloidal field components, then A^ can be determined up to a constant. The 
contours of constant A^ represent time-dependent poloidal magnetic field surfaces for any </>. Shown in Cartesian coordinates, 
beyond a few gravitational radii from even a j = I black hole, the density of lines represents the field strength in the lab frame. 
Near the horizon where the intrinsic volume of space is larger than in Minkowski space-time, the density of field lines in such a 
Cartesian plot overestimates the lab frame field strength by factors of < 2. 
For an inviscid fluid flow of magnetized plasma, the energy and angular momentum flux per unit rest-mass flux 

E = -T;/(pou'-) = -Tf/(pou') (B7) 

and 

L=T;/{pou') = T^/{pou''), (B8) 

respectively, are conserved along each flow line. For unmagnetized flows E is conserved for any stationary flow, while L is 
conserved for any axisymmetric flow. If the ideal MHD approximation (e'' = 0) holds, then the magnetic flux per unit rest-mass 
flux 

B' B^ B'f' 
$ = = = (B9) 

pQU'' pou^ poiu'^ — u'Q.f) 

is conserved along each flow line. This also implies that (<I> + /7')/po = b'' /v'' = b^ Iv^ . IBekenstein & Oronl i 1 978ft also show that 

^ = -E + ^lFL = h(ut + nFUA,), (BIO) 
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is conserved along each flow line. One can reduce E and L to forms such as 

E = -Cu,--^^b,/h = -hut + ^nFB^ (Bll) 

L = Cu^ + -^<^b^/h = hu^ + <^B^ (B12) 

where h = (po + Ug+p)/po is the gas specific enthalpy and ( = h + b^/ po is the total specific enthalpy. Notice that for an isentropic 
flow that d(po + u) /dpo = h and so 

dp/(hpo) = dh/h. (B13) 

Clearly any ratio of these conserved flow quantities is also conserved along each flow line (e.g. the energy flux per unit magnetic 
flux (£■/$)). Any axisymmetric, stationary flow solution can be written in terms of the 6 independent quantities A^, B'^, fi/r, 
$, E, and L, where the single function determines the dependent quantities B'^ and B^ . Entropy is a dependently conserved 
quantity when one writes the rest-mass density and enthalpy in terms of the entropy and another conserved quantity. In general, 
the solution is determined once the conserved flow quantities are set by the boundary conditions and the other quantities a re set 
by the equations of motion, either directly or using the Grad-Shafranov approach (for a review see, e.g., lLevinsonll2(K)l . The 
limitations of, and hence extensions to, the ideal MHD approximation are described in Meier (2004|). 

Note that for an axisymmetric stationary degenerate fluid the Boyer-Lindquist coordinate components are related to the Kerr- 
Schild coordinates by B'[BL] = B'[KS] , B^[BL] = B^[KS],B'^[BL] = B'^[KS] -BlKS](fl-2rr2f)/A, while for a general fluid 
Br[BL] = Br[KS] + (aB^)/A, Bg[BL] = Be[KS], and ^^[BL] = ^^[KS]. 

C. IDEAL MHD FLUID FORCES 

First, to investigate the spatial acceleration of the fluid along a flow line one requires a unit length space-like vector that 
satisfies V^A^.^ = since A^ is conserved along a flow line. The unit-length magnetic field satisfies these properties since 
^'''4^0, /J = 0. To study the poloidal acceleration along a field line, the toroidal magnetic field component it projected out to obtain 
the unit vector 

Bi^=NB;=N{B^'-^-^)=N{B^^-BH>;) (CI) 



where N = 1/ ^BpBpgfj,^ = 1/ yjB'BJgij with {/, j'} = {r, 6} only, and oj*" is the basis one-form. Therefore the projection of an 
acceleration along a poloidal projection of each flow line is 

flA^=fl^B; = fl,B'- + flX (C2) 

Second, to investigate the spatial collimation of the fluid, a unit length space-like vector that is perpendicular to the poloidal 
field line and perpendicular to the 9!)-direction (i.e. is required. This vector is 

C^=NcC^=Nce':^pB%l, (C3) 



where A^c = 1 / yC'^C'^ gtiv is the spatial permutation tensor Thus C = NNcB and C -NNcB''. Therefore the projection 

of an acceleration in the collimation direction is 

acoii = a^,C^'=Nc(arB'^-agB'). (C4) 

Notice that for 6 < n/lif ac„ii < 0, then the flow is collimating toward the polar axis. For > 7r/2 if 0^,,/; > 0, then the flow is 
collimating toward the polar axis. 

Third, there are many interesting frames to measure the acceleration. The acceleration away from geodesies is obtained from 
the projection of P'^" = g^^'^ + u'^u'^ on V.y7j7 = 0, giving Euler's equations for the deviation from geodesic motion 

Poha';=-P^°'p,^+rF^^ (C5) 

where "(m'^ + r''^aM'^) is the "geodesic acceleration" away from the geodesic motion (a'^ = 0). Note that F here 

is the connection coefficient and not the Lorentz factor This comoving geodesic acceleration "hides" the effect of gravity on the 
fluid. One could instead focus on the coordinate acceleration = m'^m'^, which represents the change in the 4-velocity in the 
momentarily comoving frame. From the geodesic equation of motion, 

a'^ = a'^-T''apu°'u'^, (C6) 

and so the coordinate acceleration along a flow line is 

flj^ = a^b; = a^^ -B^u^u^T^^f,, (C7) 

so the geodesic deviation and gravitational acceleration along each flow line can be studied separately. 

Notice that for a stationary flow = where the poloidal velocity = 0. For a black hole with field lines that tie the black 
hole to large radii, there must exist a region where = 0. For an outflow to reach large distances, the region wh ere u'' = 
separates those fluid elements that eventually fall into the black hole and those fluid elements that reach large distances. iLevinsor] 
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(|2005|) refers to this as the "injection surface." If particles were created only due to reaching the Goldreich- Julian charge density, 
then this must be the location where particles emerge. For an injection region with negligible angular velocity u'^ « 0, then 

fl^^ «B'T„,+B^re„ (C8) 
determines the location of the stagnation surface. In Boyer-Lindquist coordinates 

r., = ^^z(^ (C9) 

and 

2ra2cose'sine' 

^ett= ^2 ■ (LIU) 

See section|5]for a discussion of injection physics. 

C.l. Forces in Lab Frame 

The forces as written in the lab frame, rather than comoving frame, allow for simple understanding of the force dynamics. 
EQuations lC5l and lC2l implv that 

aX^=NB^{aj +aj ), (Cll) 



where for a stationary, axisymmetric flow 



is the hydrodynamic acceleration. For an isentropic flow, equation iB 1 3l implies that 



^C(MA) ^ p.,(iQg;j^ . ^ -uju'{logh)^i+ilogh)j. (C13) 



For an axisymmetric, stationary, ideal MHD fluid the electromagnetic acceleration aj^^^'' = J°'Faj/ (poh) reduces to simply 



af'=—^. (C14) 
poll 

Thu s the m agnet ic acceleration is due to the gradient of the toroidal magnetic field along a field line. Notice that from equa- 
tion or that in the limit poh that — > E/(^ilf) = L/$, a conserved flow quantity, thus BW^j = implying the 

electromagnetic field is force-free. However, notice that the acceleration aj then becomes undefined. 
Eauations lC5l and lC4l implv that 

ag,=NNcBUaf'''\af''V,, (C15) 

where e''j is the poloidal permutation tensor and af^'^'^^ is given in equations IC 1 2l or lC 1 31 This hydrodynamic collimation is due 
to the pressure acceleration in the comoving frame along a field line but directed to collimate the flow. For an axisymmetric, 

stationary, ideal MHD fluid the electromagnetic collimation acceleration (C-'J°'Faj/(poh) = C^a^^™^) reduces to 

^,jEM, ^ NNc (,«^(b.)2^[5^] + ,^,Vb*B,,,) , (C16) 
Poh 

where f[B„] = {Ba,b -^Fij aBg).!,^ and (B'')^ = {B'')^ + (B^)^ and e„'' is the poloidal permutation tensor The last term on the right 
hand side of equation IC 1 6l represents the "hoop-stress" that leads to collimation for nonrelativistic winds. The first two terms on 
the right hand side correspond, respectively, to the forces due to poloidal magnetic stresses and the electric field (£') gradients. 
The latter can collimate relativistic outflows. 

D. COMPTON SCATTERING 

In the lab frame, seed photons are Compton upscattered if the energy of the photon Eseed <C TgineC^ for an electron Lorentz Tg. 
The upscattering continues until the lab frame photon energy exceeds electron energy. In the lab frame, each scatter gives the 
photon a new energy Escat ~ ^^^e^seed if E^eed < meC^/^e and Escai ~ TenigC^ otherwise. Two limiting scenarios are if the photon 
crosses the jet or if the photon travels parallel and within the jet. 

The optical depth to Compton scattering for a photon in the rest frame of the jet electrons is 



PO,e-e+ 



ajdl', (Dl) 



which in the lab frame gives 



(D2) 
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where /3 = v/c, = (1 - (3^) 'Z^, and 9 = 11/2 corresponds to perpendicular interactions and 9 = to parallel. Across the jet 



while for along the jet 



where /3 = (^f^) . For T, > 1, r,(l -/?) w (2r,)- 



r± = / ( I aTT,rd9, (D3) 



PO,e-e 



aTr,{l-p)dr, (D4) 



where oj is th e Thomson scattering cross-section and ro is some fiducial starting radius in the jet ("see lRvbicki & LightmanI 19791 
lLongaiiil992l) . 

Similar calculations can be used to estimate the forward or backwards pair-production or pair-annihilation optical depths. In 
the lab frame, for a photon gas moving with F with one photon energy E = TE' and another > (FmgC^)^//?, then the photons 
annihilate with a cross section at £ < O.SMeV of a ~ ctt- and the Klein-Nishina corrections for photons with E > O.SMeV modify 
the cross section such that a oc E~\ For a spectrum (number per unit time per unit area per unit energy) fE'"", then for a = 2, 
the average cross section is cr w 0.06cr7-. In this case the relevant proper density is when for a photon beam, where n-y is the 
number density of (typically fewer) high energy photons. See also Lithwick & Sari (2001). 

For a beam of electrons with velocity j3, a / P for nonrelativistic electrons and for F ^ 1 

a«^(log2F-l). (D6) 



The pair annihilation rate is 

-1 ^ / \( 

which can be compared to some dynamical time to determine if pair annihilation is important. 
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